4 


VOLUME 87 NO. IR2 

«6 (ep AMERICAN SocteTY 

NY OF CIVIL ENGINE 


—— 


reach ae is given asa with paper. "Those who are pede 
l 


of submit material will expedite the review and f on" a rocedures by com 


f approximate y 0 words par pi 
en of conclusions ‘must end it. 


. ‘The ribbon coy copy co} 


of all illustrations, tables, etc., must be 
The author's full name, Society membership grade, rade, and 


stating present em the first page of the pape: 


Mathes ics rece he copy that is submitted. Because 


is, it i be drawn carefully, ey that special symbols be 
¥ roperly identified. The letter symbols used should be defined where 
appear, in illustrations or the text, 


{an origina 
of 8 by 11- 
on on must be made in the text for "each _ 


. Ilbustratiors bé drawn in black in 


per. Because illustrations will be with a width 
nd 414-inches, the lettering must be large enough to be legible at this width, — 
hoto: raphs s should be submitted | as glossy prints. and 
- made within the text for each illustration, 
The desirable average a paper is about 12,00 000 and 
absolute maximum is 18,000 words. As an approximation, each ac page of typed 
text, table, or illustration is the equivalent of 300 words. 


this be made on condition that th 
of author, name of publication, page reference, and date of 
_ Society are given. The Society is not responsible for “ay catemens made or ‘opinion 
in its publications, 


‘ion Publication office is at 2500 South State Street, Ann 


"8 . Editorial and General Offices are at 33 West 39 Street, New York: 18, | 


$4.00 of a member's dues are applied as a subscription to this Journal “Second: 


class postage paid at Ann Michigan. 


BASIC REQUIREMENTS 
f t rs should be submitted to the 
Manager of Technical Publications, ASCE, Authors sho 
— 
eference 
— 
— 
| 
— | 
“Bibliography,” which would encompass the su 
— 
indexes for previous years are also 


Sournal of the 
-IRRIGAT ION, AND DRAINAGE DIVISION 


> IRRIGATION AND DRAINAGE DIVISION 

EXECUTIVE COMMITTEE 
Kenneth Q. ‘Volk, Chairman; N. AL | Christensen, Vice Chairman; 
? Paul H. Berg; Herbert E. Prater; William W. Donnan, Secretary = 


Merrill Barber, Board Contact Member 
OMMITTEE ON PUBLICATIONS 


Dean . Peterson, Jr., Chairman; Charles L. Barker; Harry F. 'F. Blaney, S Sr.; 
_ Eldred R Harrington; Sol D Sol D. Resnick; Daryl B. Simons; ‘Calvin C. ‘Warnick 


June, 1961 


& 


Field Measurement of Conductivity 
W. Donnan and V. S. Aronovici . 
oO 


Outlook. for Economic Use « of Fresh Water from the Sea 
Unsteady Flow of Ground Water into Drain Tile 
‘Water Table Fluctuations” Induced by Irrigation 


Copyright 1961 by the American Society of Civil Engineers. 
_ Note.—Part 2 of this Journal is the 1961 -21 Newsletter of the Irrigation - 


7 three issues of this Journal are dated 1960, 


Proceedings of the American society Of Civil Engineers 

— 

| 

7 

&g 


‘World in Wi Water at 
Charles WwW. Thomas. 


December, 1960. Discussion closed. 


(June, 1960. Prior ¢ discussion: 


Los Angeles Water Supply and Irrigation, by Samue Samuel B. 


Morris. (December, 1960. Prior discussion: None. 
Discussion closed.) 


Method for Estimating Use of. Water f for 
_ Agriculture, by Wendell C, Munson, (December, 1960. P 


_ discussion: None. Discussion closed. ees 
by George H. Hargreaves” 


“Irrigation Systems of the Tigris and ‘Euphrates Valleys, 


‘Stanley Ss. Butler. (December, 1960. Prior discussion: None. 


by Vahe J. Sevian and Kamil Taj- 


L* — 
— 
tm 
— 


Journal of the 
IRRIGATION AND D DRAINAGE 
ns Proceedings of the , American Society of Civil E Engineers _ 


FIELD MEASURE han OF HYDRAULIC ConpucTIVITY® 


‘By’ William W. 1 fF, ASCE, and V. V. 8. Aronovici 


4 ‘One of the major factors ‘connected with the éenten of drainage systems is 
the measurement of hydraulic conductivity of the soil strata to be drained, This 
paper describes the work carried on to develop a device for this measurement. 
-— Laboratory studies resultedinthe adoption of a small, brass, screen-type well 
_ point which could be inserted in the soil below the water table. The rate at 
which water pumped out of the well point related directly to the | -conduc- 
Field-type well points were developed, using ordinary hardware pipe > 
_ fitting. The field well points were then n tested in th the laboratory | and in the field J 7s 
on _ This report describes how to make the field well points and wine outlines a 
- technique for for field installation and operation of the device. In addition, graphs 


se Note.—Discussion open until November 1, 1961. To extend the closing date one 
month, a written request must be filed with the Ensoutive Secretary, ASCE. This paper 
‘ts part of the copyrighted Journal of the Irrigation and ee 0 Division, Proceedings 
of the American Society of Civil Engineers, ‘Vol. 87, No. IR 2, June, 1961, 
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INTRODUCTION 
The Soil Conservation Service (scs) has been _—e upon to design drain- _ 
age systems as a part of their overall program of technical assistance to Soil . 
_ Conservation Districts. One of the major | factors connected with the design of 
system is the measurement of hydraulic conductivity of the soil 


— 


water | will move to the drain in order to space the drains intelligently an and to 
_ determine the speed with which the drains will dewater the soil profile. The 
_ SCS has listed this problem as one ofits high priority research needs. There-| 
- fore, the objectives of this project were to develop a relatively simple device — 
=: in measuring hydraulic conductivity rates of sand strata beneath the 
water table, to test this device in the field under operational conditions, and to 
set forth a “eee and neuatins which could be used in the field by f field 
EXISTING METHODS AND TECHNIQUES — 
Various methods, devices, and techniques have been proposed for ne 
_ameasurement of hydraulic conductivity of soil strata in place. The most sim-_ 
_ “ple method is the single auger- — method proposed originally by E. . Diserens3, 
‘Sage, 6. B. Hooghoudt 4 improved on earlier developments and = 
= F. Ernst°, using the relaxation method, developed an empirical formula 
which is claimed to be Field were developed by C. H. M. 
van Bavel and Don Kirkham®, OW. C. Visser’ later simplified and developed a 
a field method which | is now widely used in Europe. A detailed are of this” 
technique has been completed by W. F. J. Van Beers®, now = 
Field methods other than the simple auger-hole method —_ —o the 
use ofa piezometer have been proposed by Kirkham9, Frevert and Kirkham10, — 


a 3 “Beitrag zur Bestimmung der Durchiassigkeit des Bodens in naturlicher Boden- 


lagerung, ” "by E. Diserens, Schweiz, Landw. -Monatshefte, 12, (1935, Pp- 188-198, 

tot de von natuurkundige grootheden van den Grond 
B. Hooghoudt, 4 Versl. Landb. ond., Vol. 42(13)B, 1936, pp. 449-541.” 

_ 5 “En nieuwe Formule voor de Berekening van de Doorlaatfactor met de Boogaten- 
= methode,” by L. F. Ernst, Rap. Landbouwproefsta en Bodemkundig, Inst. T.N.O., Gron- _ 
6 “Field Measurement of Soil Permeability using using Auger r Holes,” 
- Bavel and Don Kirkham, Proceedings, Soil Science Soc. of America, Vol. 13, 1949, pp. 
“Title Drainage in the Netherlands,” by | W.C. Visser, Netherlands Journal of Agri- 

cultural Science, Vol. 2, 1954, pp. 69- 
"Auger Hole Method,” by W. F. J. Van Beers, Bulletin No. Internat. Inst. for 
Land Reclamation and Improvement, Wageningen, Netherlands, 19 58. 
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Table,” ‘by Don Kirkham, , Proceedings, Soil ‘Science Soc. of ‘America, Vol. 10, 1946, 


10 “A Field Method for Measuring the Permeability of Soil bikew a Water Table,” oy 


4 R. K. Frevert and Don Kirkham, Proceedings, 8, Highway | Research Ba., Vol. 28, | 1948, ‘Pp. 
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and Kirkham!1, ‘and Reeves and Kirkham12, In. addition, E. Chitds13 


proposed the use of a two- -well method, where water is pumped out of one well = 
= into an adjacent well. Subsequently, this technique was a upon by 5 


utilizing three- -well and four-weli methods. 


wan te all of these techniques, while valid, are of necessity rather cumber- 


a 


All of the known techniques for measuring hydraulic | conductivity of the soil 
; _ strata in the field have one common weakness. This is that the “A” function, 

or area and geometric shape of flow into the auger hole, piezometer, well, or 
other device is difficult to determine. Furthermore, once itis determined, this 
factor ‘is difficult to maintain under most field conditions. It is a well-known © 
- fact that flow of water into an open cavity in the soil creates sloughing which, 
in turn, alters or changes the size and shape of the cavity. The present work 
: to date (1961) has been aimed at developing an instrument with a fixed or te 
determined “A” function, thus this troublesome variable from t the 
The initial objective also was ~ develop a relatively economical unit | 
could be readily ecient or made without resorting to elaborate laboratory | 


device | in . Fig. 1 is a sintered glass gas diffuser; the center device | isa ( ceramic 
fish tank airstone. The technique used was to insert these devices in soil to _ 
tested andto apply a low head suction to the outlet tube. It was quickly apparent 
the airstone was not applicable because of the tendency to clog. Se oe 
_Aseries of tests were conductedusing filter paper wrapped around the iis 
; stones, but still the results were not satisfactory . The airstones tended 7 


age ¢ andovera periodof several days the stones themselves would ceevenss in 


permeability even in clear water. "ve 


_ “TI “Piezometer Method for Measuring Permeability of the Soil in Situ Below a Water A 
Table,” by J. N. Luthin and Don Kirkham, Soil Science, Vol. 68, 1949, pp. 349- 358. 
12 “Soil Anisotropy and Some Field Methods for Measuring Permeability, ” by R. C.. 
Reeve and Don Kirkham, Transactions, Amer. Geophysical Union, Vol. 32, , 1951, 

“Measurement of Hydraulic Permeability of Saturated Soil in Situ. 1. Series A, 


7 of a Proposed Method,” by E. C. Childs, Proceedings, Royal Soc. of Lenten Series A 
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. | The next step was to construct some type of a brass screen well point which J. i 
_ would havea fixed or constant “A” function and which could be used in place of = 
an ordinary piezometer or auger hole. Several experimental well points were - aa 
: tested. These included double core typed wherein the inner core could be — 
| 
= = 
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_—-WELL-POINT DEV ICES USED IN HYDRAULIC 


— 
2,—NON-C LOGGING BRASS SCREEN WELLPOINT 
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"wrapped with a final design, and the one with which m most of 
_ the subsequent laboratory t testing was carried out, is shown as the a “hand 
device in Fig. Fig. 2 shows this ‘point, unassembled. 


“forated brass tube with an outer cage assemble and demountable brass 
The outer cage is covered with 40-mesh brass screen. The test point has a 
screen length of 3.5 cm (1.375 in.), an outside diameter at the screen of 1.85 cm 
(0.728 in.), and an overall intake area of 20.29 sq cm (3.146 sq in. ) a ve 

7 _ The “A” function for this device was originally c ‘computed after Luthin and a 
Kirkham!1, They evaluated the “A” function by means of a three- dimensional 
; electrical analogue of -” ground-water problem. . However, this work was done a 


Cavity 4 in. long (Open end-Luthin 


and Kirkham) 
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Diameter of in inches 


_GRAPHS FOR COMPUTING FUNCTION. 


“usinga cavity with an open-bottom end whereas the test has a closed 
end. Subsequent investigation by N. A. Kadirl4 revealed the existence of the _ 
desired “A” function. Kadir developed data for computing the “A” function val- : 
-ues for sealed bottom cavities. . The Kadir cavity more nearly approaches | the 
geometric boundaries of the Pomona well point. Thus, the Kadir curve has been 


used to compute function based on data by Luthin and 


“Measurement of Permeability of Saturated Soils Below the Water Table, by ‘N. 
Kadir, Ph. D. Di ssertation, Utah State Univ., » Logan, Utah, 1955. 
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une, 

‘An example of how to compute the “A” functio n is given as follows: The : | 

- ratio of length to diameter of the brass screen device is computed to be 1.0 to E : 
0.529. The graph in Fig. 3 for computing the value of the “A” function is drawn — » 7 

for a 4-in. length cavity. Thus, a diameter to maintain geomet- 

: ric similarity would be 4.0 x 0.529 in. = 2.15 in. diameter. Referring to Fig. — 
3, the “A” function value after Kadir would be interpolated as 15.3 in. fora — 
4. -in. length “ee The value, then, for the brass screen device of 1.375-in. 7 
length would be (1.375/4) 15.3 = 5.25 in. This value has been used for all com- 

: _ putations involving the brass screen well point. The value is assumed to be 

. validfor any measurement of hydraulic conductivity made at a reasonable depth = 

below the water lot 5 in. or more below, and for a reasonable height — 


s screen ~~ point v were conducted in a tank 


45 cm high, and 95 cm long. The front face was fitted with : ail-cm thick, clear 
_ plastic plate. Water supply reservoirs were provided along the full width and | 
height of each and with access of water into the entire soil profile through 
, _--« Sereened access holes. Provision was also made for input and withdrawal of © 
water from the bottom of the tank through a 3/4-in. _ diameter perforated pipe = 


ne along the floor of the tank and outletted ad through one end. sesified as 


was installed in 1 the approximate center of the tank at a Po ‘of 20 cm below | 

the surface of the sand in the tank. Installation was made by pushing shing the well 
point vertically into the saturated sand. 

a Suction head v was: introduced to the well point by m means of a siphon which 

- could be adjusted to provide various head differentials for various test runs. 
Water levels ‘is in the tank ¢ and suction head on the well point were both m measured © 

Si from the same meter stick affixed to the outside front face of the tank. ‘The 
_ suction head was reflected in a glass manometer rae along the side 


of tl the meter stick. 


ih making a a test run run of the well point, the sand in ‘the tank was comatetels 
saturated by the introduction of water from the floor intake tube. The initial 
test runs were not made until the tank has been saturated and dewatered to field © 
| Sli through a number of cycles. This allowed for a reasonable degree of 
2 settling and repose of the sand prior to actual well-point tests. = ~~ 
_ The water in the tank | was then brought to : a ponded depth of approximately “- 
1 cm above th 
ing a test run, water was introduced from the end supply tanks fed through a 
Mariotte arrangement. The pond surface was heldat a fixed level | level by means of 
a strong ‘suction outlet suspended at the de desired pond level. i fom abel 


— 


: 
| 
| The laboratory tests of the bras: 
— 7 
Two manometers were provided by embedding ceramic airstones at approxi- 
 &§ g mately the third points within the tank. These airstones were fitted withtubes §f 
RUNS 
— 


removed from the siphon to flow out of the tank 

- through the siphon with the well point acting as the i intake source. ‘The volume —_ 
= of flow through the well point then became a function of suction head as provided — 
by the siphon tube. Since the tank was mounted on top of a 4-ft-high laboratory 
= the siphon arrangement provided for a convenient range of of suction head 


0.1 cm to about 80 cm. 4 


was ina beaker r reading to 5cucm. “No tempera-— 
ture corrections were made since no attempt was made to provide « continuity 


between tests fro hour to hour or day today. a 


“TABLE 1, _-HYDRAULIC CONDUCTIVITY OF HANDFORD FINE SAND AS MEASURED = 
BY LABORATORY PACKED 3-IN, CORES 
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The hydraulic conductivity of the Hanford fine sand in the tank was checked _ 


_by the use of laboratory- y-packed c cores. These were the conventional 3- in. di-— 
ameter cores 3 in. long. Table 1 shows data from a series of permeameter — 
— tests of soil cores. The apparent Fccsaaes conductivity. of this sandis between © 
45 20 in. per 


| 

| 

No, 

204 500 500 31.48 [9.78 18.9 
204 | 500 | 500 17.8 

252 925 500 175 2,82 | 22,9 
10 500 | 490 15.79 | 8.19 | 22,1 | 

252 | 50 450 8,38 | 
252 | 500 | 500 26.55 | 9.05 | 20.4 
252 | 500 500 8.60 21.5 | 
250 | 505 | 255 1048 | 
201 | 500 26.65 | 1219 | 152 
| 
the earlie point in the Hanford sand indi- 
cated that cons Thanifested in the apparatus used 
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to provide and record suction to the well point. To circumvent this and to make 
head loss measurement, anauxiliary manometer was installed inside the well 


point itself and outletted via suitable tubes to the front face of the tank. This 


_ provided for the ‘simultaneous recording of both the siphon tube ‘manometer 


a and the apparent suction head, as manifested by the water level in the well- “i 
: point manometer. Table 2is a typicalexample | of test runs with the well ell point. — 


TABLE 2. —OBSERVED DATA BRASS WELL- _POINT TESTS 


Well. 
point» 
| mano-— 
meter, meter, — 
inches | inches 
per 
hour 
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Well-point manometer 


should also be noted that when the applied head was used in successive tests 
on a descending scale of magnitude (trials 16 to 21), the comparison between 

successive conductivities was more favorable. In other words, if the well point — 
stressed, trials with differential heads did not vary 


| 


— | 
“hg, “hg, | hy -hg,| T, in] Q, in 
— em? | cm om | cm | utes q mi 
.2(] 10.70 | 13,81 
— 
5 10.38] 13.85 ff 
— 
t+ 
— | 500 | 345.0 | 9.40 | 13.50 
— 480 | 240.0 | 9.40] 13.100 
— 400 | 200.0 | 9.50 12.720 
415 | 138.2 | 9.57 | 12.43 
265 | 88.4 | 9.94] 12.41 
132 | 44.0 | 9.55 | 12.39 
raulic conductivityas determined 
Cores 4nd the conductivity tas Pom the wo Doint toct 


_— phenomenon prompted a series of tests wherein a given suction head 
=e appliedand allowed to run overan extended period of several hours. Peri-_ 
odie readings were made to record drift or fluctuations in quantity of flow. It 
was found that after about 30 min the flow came into equilibrium. Subsequent 
_ tests w were made wherein the well point was stressedat some higher given suc- a 
tion head and then the head was to th the de: desired amount. The 


in. standard galvanized nipple with 20 2 ig in. 
perforations 


2 in. nipple - 


4.—POMO} NA WELL POINT 


flow equilibrium. This suggested a technique which 
was found to be advantageous when the device was tested in the field. phe 
_ Again referring to Table 2, it is apparent that excessive suction heads cre- 
ated excessive head losses. These occurred not only in the siphon manometer | 
ae but also to a degree in the well point manometer observations. 
7 > This circumstance, while nota factor in actual field applications, does suggest 
4 ) _ the use of low head differentials infield use. Low head differentials would tend 
_ to minimize the head losses which will be a | characteristic of aan — of we 


1 

t= ==2 |, in. Standard galvanized pipe- tap thread 

| 

4 

4 

q 


4 through the brass screen of. well point. At present, | the 
for field use call for for only a 3-in. or 4-in. head to be to the a 


FIELD | WELL POINTS LL POINTS 
he primary of this project was was to instrument for 


Cs technicians and others. Thus an attempt | was made to construct a — : 


4 


Hydraulic 
-incoretests | 150-lbdrumtests | Well pointtests 
po 


umber per hour | per hour* i | perhour 


718 | 59 | 7/20 
7/18 ry 7/21 
7/20 | | 
7/21 | 7/25 

7/26 


3 


4 Value is arithmetic average of 10 individual trials as shown in Table 4, 


well points consist of ordinary galvanized pipe nipples, couplings, and plugs © 
and are constructed as follows: Cut a 6-ft length of 3/8-in. galvanized pipe 
tap thread the inside | ends with a in. ‘standard pipe ‘For the inner 


| 
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oH 
TAB 
type 
(ly) 
— Antonio | 15 7/18 | 13 7/21 18 
Average | 152 | | 152 | 165 — 
7/20 if 24 | wa 
+ 
— 
| 


: pe solder it 1 into a cylinder 3 in. plato screen slips loosely over the 2 


Be For the tip-end con component, use a standard 1/4-: in. . galvanized pipe Jf 


_ with a standard 1/4-in. pipe plug threaded on one end, _ The end peice threads 
_ onto the end of the nipple and should be threaded up snug against the cylindri- 
4 cal brass screen. The blunt tip and plug can be shaped to a conical point. ‘The — 
TABLE 


PP ‘pumped, in | time, in cu cm 
permin jin inches 


— 109.0 


109.5 


2 
5 


109.7 


4q x 2 In tests carried out under a low vacuum, the pump used to remove water from in- aad 
side the well point was set at a speed which just barely | _ jun of the inflow | 


“side the well was set at a considerable of air along © 


field- -type point ‘costs about $1.25 25 and fits through the end of a 


standard Viehmeyer King tube. See Fig. 4 for a —- of this device. uy im 


‘TESTS OF FIELD WELL POINTS pe 
The field- -type well point hi _ been tested extensively in the laboratory. The 


echnique used was to pack approximately 150 lb of sand material into 4- ft- -high — 
steel drums. The drums were provided with water inflow and outflow pipes. 
and were fitted with manometers sothat they could be used as a permeameter. . 
‘Thus the sand material in the | drums could be for hydraulic conductivity 


m 4 HYDRAULIC CONDUCTIVITY | 

2en 

> 

ame 10 200 1.84 | 1085 | 13,5 

— 

0 


aa into i a with the well point about 12 in. from the surface. A : suc- 
tion tube was then lowered inside the well point shaft. The end of this tube was — 
_ placed 3 in. below the ponded water surface. When suction was applied, a 3- in. 
hydraulic head was exerted onthe well point. Suction was provided with tha suit- 


= 


” function =7.5 in. 
Pomona 


K, in inches per hour 
‘FIG. _RATIO BETWEEN q AND K 


able vacuum quantity flow was measured. From these data a the hydrau-— 
lic conductivity could be computedanda comparison made with the conductivity = 
as measured by the 150 lb, 4 ft high drum, permeameters. A third compari-_ + 
son was obtained by preparing small core samples of the sand material and 
testing them for hydraulic conductivity in the 3-in. permeameters. . 


Mad 
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HYDRAULIC CONDUCTIVITY | 
a Four different types of sands, having a | range of hydraulic conductivity from L- 
_ approximately 15 in. per hr up to 75 in. per hr, were used inthese tests. Table — 
3 shows the results. Table 4 is a series of observations ‘made in the labora- 
. tory using the well point in a Hanford fine sand. ee iia 
as To compute hydraulic conductivity, an adaptation of the well-known Darcy 
Law of the flow of water through a saturated porous media has been used. . 


x represents the hydraulic gradient; and A is the area of flow. — 
7 ‘The adaptation as developed by Kirkham is as follows: 


er 


4 


shows the between the of water pumped for vari-_ 
ous suction heads applied and the corresponding hydraulic gradients. These 
curves 's have been 1 computed for use with the field- ~type Pomona well Point, hav- = 


well point in in the fieldc can be made by driving the soil tube 
to the approximate depth where a a measurement of hydraulic conductivity is to 
a made. Then, by inserting the well point in the empty soil tube, both devices 
= lowered into the hole. The well point t is 1 then pushed on down 6 in. or 8 in. 7 
=a the soil tube into the stratum where a measurement is to be »made. — —_ 
Mi water table is allowed to come intoequilibrium and is measured. Then 
a small-diameter suction t tube is ] lowered inside the well point pipe to a point 
3 in. below the water table. By ‘pumping o or providing a suction to this tube, 
a 3-in. hydraulic head is applied to the well point. Pumping can be done either 
by an attachment to the vacuum on a pickup tr truck or byanordinary F Ford fuel 
_ AFord fuel pump can be fixed toa 2-in. by 6- -in. board and can be saaiiaiaiah 
by hand. This device produces a 9-ft to 10-ft vacuum which will suffice for 
most field observations. The quantity | of water pumped ina given period of 
‘time is determined. Usually it is best to pump for about 5 min to 10 min and 
then take several short-term readings of rate of flow. Using the curves in 
5, the the hydraulic conductivity can be by diret ratio. 
we CONCLUSIONS | 
OA small, screened well point has been developed for use in dite measur-_ 
‘ments of the hydraulic conductivity of sand strata below a water | table. ~ 
- device has been tested in the laboratory and in the field with good results. 
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' _ OUTLOOK FOR ECONOMIC USE OF - FRESH WATER 
16 


al The need and ew for fresh water from the sea or from inland saline 
waters are mainly dependent on the cost of producing such water, storing, and — 
4 transmitting it to places of use in competition with « other ‘sources of fresh 
waters 
Foreseeable costs of “conversion far exceed of major 
_ domestic , industrial, and irrigation water supplies. Typical costs are pro 
4 sented for municipal, industrial, and irrigation supplies and for current a 
foreseeable costs of fresh water conversion. Comparison is made of water 
- availability and requirements in the humid eastern United: States a and the es 
- Conversion costs are too high for major sources of supply but are ne 
s for isolated plants and, possibly in the future, for larger general use on 


_ the western Gulf coast of Texas, which will ultimately be short of water. — 


ee -Note.—Discussion open until November 1, 1961. To extend the closing date one 
month, a written request must be filed with the Executive Secretary, ASCE. This paper 
- is part of the copyrighted Journal of the Irrigation and Drainage Division, Proceedings 
of the American Society of Civil Engineers, Vol. 87, No. IR 2, June, me pre 
_ @ Presented at the June 22, 1960 ASCE Convention in Reno. Nev. — ify 
1 Cons. Engr., Los Angeles, Calif. 
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version: of inland saline waters” is further faced with | expensive ve disposal of 
P wastes, a subject that should be under, study in parallel with desalinization © 


waters. are mainly dependent on the cost of producing such water, storing, 
and transmitting it to places of use in competition with other sources of 
- fresh water supply. In fact, even if nofresh water is available as a competing 7 
- factor, high cost of such fresh water f from the sea or from saline waters in — 
any given area not justify an industry, community, or agriculture in 
_ locating in such area. Conversely, other conditions may dictate the use of a 
given site for mining, industrial, or military purposes almost irrespective of 
cost of water necessary to the establishment and its operating per- 
sonnel. Quite obviously, these latter conditions affect special limited sites and _ 
_ developments and are not in in the nature of general urban 1 development or ir- ; 
a "The purpose of this paper will be to consider the costs of normal water — 
‘Supplies | in comparison to present and anticipated future costs of producing © 
fresh water from the sea and from saline waters; then, to relate this cost 


> -qgeeal to 0 possible needs within the continental United States, except 
oF 


an WATER NEEDS ARE DEPENDENT ON WATER costs 
: - The need for water is not independent of cost, especially for irrigation 1 use, 
If adequate high quality water costs only $1 per acre-ft, all semiarid lands of — 
proper elevation , topography, and soil type willbe irrigated when markets for 
their crops develop . However, if irrigation water costs as much as $20 per > 
 acre-ft (6.14 cents per 1000 gal), then only favored lands and locations can 
‘sufficiently high-value crops, such as cotton, citrus, fruit, nut, 
truck gardens, to pay out. In general, municipal water-users have been able | 
to pay whatever the water costs. The cost and character of water use, however, — 


become an important factor in the location of an industry requiring large iad 


of water in relation to the market value of its products, we a 
are generally of the order of $1 to $10 per acre-ft (0. 307 - e 3.07 cents per 
1000 gal) and may go as high as $40) per. acre-ft (12.28 cents | per 1000 gal) for 
“limited areas Ofhigh-value crops, 
on While retail costs of municipal water delivered through the customer’s 
meter generally vary from 10 to 30 cents per 1000 gal ($32.60 - $97.80 per 
_acre-ft) the water production costs are generally 3.07 to 6.14 cents per 1000 
al ($10 - $20 per acre- mag for local ‘supplies and 6. 141 to 21. 49 cents per 1000 


Senate Select Committee on National Water Resources Reports on Water 


- Costs.—It —It will be be helpful to examine the prices that are being paid for water — . 


veyed hundreds of miles through large aqueducts. 


&g 

| 

: 
— 
i 
it 
| 
— 
| 


FRESH WATER 


reported to the Select Committee on National Water Resources, United 


“With ‘respect 1 to the cost | of water for industrial use, there is a a 
variation. The major water- -using establishments ‘usually self-supply 
their = requirements, attempting to hold the cost within the range of 5 a 
- 10 cents per 1,000 gal ($16.30 to $32.60 per acre-ft). On the other hand, 
smaller volume -_water-using establishments which | average only 
from 4% to 16% as much water, usually purchase from a public water 
utility | at costs ranging from 12 to 20 cents per 1, 000 gal ($39.12 lll 
— $65.20 per acre- -ft) for larger volumes» and from 20 to 28 cents per 
i 000 gal ($65.20 to $91.28 per acre-ft) for the extremely small volume 
a users. No comprehensive assembly has been made of these costs, = = 


form to usage, 
These prices obviously are waterworks sale prices and not water 


Bore noy for “smaller volume water-using establishments,” $13, 04 to $21.73; 
and for “extremely small volume users, ” $21.73 to $30.43 per acre-ft, These 
a figures indicate the quoted cost to major water-using establishments is out of 

line for their own self-supplies which must be materially lower than prices 7 

paid to public utilities. Experience has shown that major v users generally ob- - 

_ tain their fresh water supplies at at $5 to $10 peracre-ft, 4 4 


a actual consumption or or depletion « of water asa resource is indi- 

cated in the 1954 Census of Manufacturers inthe difference between the 
intake and the discharge of the manufacturing industries. In aggregate, 


depletion is less than 10% of the total 


reports of the Federal Power Commission (FPC). It is unfortunate that there 
is an absence of such accurate reporting of water production, water sales, and | 


& hour produced and sold are available for the entire industry in the aie 
By ay use of water. It is recognized that such data like hydrologic data 


z more difficult to procure, and involve wr estimates from incom- 


4 into the false belief that there is great impending water shortage if fresh water 
not obtained from the sea and saline waters. 


Present and Anticipated Future Cost.— Turning» now to the present (1961) 
“and anticipated future cost of producing fresh water from the sea: 
_ Senator Kerr’s letter of July 13, 1959, to the Secretary of the Interior, to 
oe attention of the Director of the Office of Saline Water? | contains, in part, ‘ 


a a The report should contain a brief description of desaliniza 


cesses now in use, both here and abroad, discussing engineering p princi- 


as so-called This misleads the public and their 


aed 2 Committee Print No. 8, April 1960, under the title, “Future Water Requirements oO : 


Principal Water- Using 16. 
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ples, co: costs pi per thousand gallons and per acre-foot, and environmental | 
_ conditions, and actual areas within the United States suited ng the appli- | 


cation of these processes. - 


The Senator Kerr letter and asks for similar information on 


rocesses in the ilot lant st stage, and the we ee 
pr in the pilot plant stage, andthe 


“best to desalinization will be practiced by 
1980, giving geographical areas and amounts of water that can be pro- 


2s The reply is not quite as detailed as the ‘request. _ Following is the reply by | 
"the Assistant Secretary of the Interior on October 19, 1959, _ writer has 
4 i ivalent cost it feet): 

nserted in each ach case | the equivalent cost in acre- 


yous 


«Just prior to World War the cost of 


one or more > of the plants Public Law 85- - 833 authorized the iret 
a to construct and operate, we expect to. convert saline water to fresh for — 
a less than $1 per thousand gal $326 per acre-ft)--- New © fresh for 
*y 3 still in the laboratory stage and not advanced far enough to be considered - 
for demonstration plant programing, give promise of lowering the cost 
: of conversion below 50 cents per thousand gal ($163 per acre-ft).” 
 #For years the citizens of Coalinga, California, like many other com- — 
hs munities in this country have had to haul in their drinking water supplies a 
In Coalinga, this water was obtained at a cost of $7 per thousand gal 
— ($2282 per acre- ft). Earlier this year this community cut its water bill 7 
to $1.45 per thousand gal ($473 per acre-ft) and made U. S. history by 


7 Pre the first to get its drinking water supply from brackish well 


Elsewhere from this’ report2 there is extracted some additional cost data 


ss “The cost of fresh water produced by the most efficient existing sea-_ 
mate conversion plants is now approximately $1.75 per thousand gal 
7 ($570 per acre-ft). It is anticipated that the first two sea-water distilla- 7 
os tion demonstration plants will produce fresh water from the sea for $1 _ 


or less per thousand gal ($326 per acre-ft).” 
_ Estimated costs of producing fresh water from inland saline waters were 
_ In this table, under column 2, salts are increased from 250 ppm to 500 | ppm 
for “Projected Intermediate Size Plant,” thus defeating comparability. All 
costs are higher, however, than the irrigator can afford to pay. 
From the preceding data on water costs, it is apparent that even projected - oy 
costs of obtaining fresh water from the sea or from inland saline waters are 
too high for agriculture. Such methods do not afford a hopeful substitute for 
water obtained by normal methods of — storage and conveyance in in long 


LITTLE CONSUMPTIVE USE OF 


In general, the areas of the United States east of the Mississippi that have | 
= 40 in. or more of rainfall suffer little saat use of ° water, except along 


| 

| 

& 
| 

— 

f : 

4 
> 
| 
— 
4 
— 
| 
= 


(ead ‘000'r +) 
og" 


LINN Xd TVALOV 


quouwr 


 SUALVM ANITVS CNV'INI YALVM ONIONGOUd AO SLSOO aa.LVWILsa—" 


ee 
q 
: 
3 ! 
4 a 
q 
|) —= 


oast lines where sewers carry domestic and industrial wastes to the sea. 7 
The major ‘problem is. that of the treatment of these inland wastes and the — 
treatment of raw water before re-use, as our lakes and rivers afford both the 
supply of water and the conveyance of wastes. Such far below wthe 


7 foreseeable costs of pr producing fresh water from the sea, ; 
ant 
SEMIARID 


B For the reasons _ in the preceding section this paper r will confine itself 


and semiarid areas of the United States. It contains, also, along the northern ~_ 
Pacific coast ranges, the highest rainfall (more than 100 in.) and much of the - 
high mountain areas of the Cascade, Sierra Nevada, and Rocky Mountain ranges © 
where rain is generally in excess of all vegetation requirements. Such sal 
are the principal sources of the streams and rivers of the West. 
aes With the exception of the cities and areas of the Pacific northwest, all ee 
areas in the West, including nearly all cities , metropolitan areas Of major _ 
population and industry, and irrigated agriculture lie principally ir in the regions 
mean | annual r; rainfall is low, _ varying from 5 in. to 20 
tal United States, 
éiiethinan Alaska, but average only 27% of the forty-eight states’ water runoff. 

They include, however, in the Pacific coast states, Arizona, and Nevada, the 
most rapid growth in population and industry inthe nation, with the ‘single ex- - 
ception of Florida. ‘Nearly all irrigated agriculture is inthe West. 

It is here, then, in these seventeen western states, that water supply and 
waste disposal problems are of truly major proportions. Here are areas of 
least rainfall and lowest runoff. The locations of irrigable lands, population | 
and industry in relation to places of large rainfall and stream flow are most J 

_ displaced. _ There are vast areas (of closed basins such as the Imperial Valley 3 
ae of California and the Great Basin centering in Nevada and Utah. And yet, here . 
are the nation’s largest dams, reservoirs, hydroelectric plants, longest aque- — 
ducts and greatest per capita consumptive v use | of water, The unfamiliar will 


say that here the need for fresh” water from salt and saline. waters. 


older, more e populated a1 areas of the United States and of the old world have not | 
yet been made. The population has the vitality, spirit, and optimism of youth. 
Where else would one find such a bold proposal as the California Water Re-_ 
sources Development Bond Act for $1,750,000,000 authorized a the last session 


their vote at the general election, November 8, 1960. <egeineene el 
‘This project will provide 4,000,000 acre-ft per annum, three billion five je 
hundred million gallons of water per day, to be stored and transported as far @ 
as 600 miles through a great multiple-purpose water project. Flood control, 
_ salinity control, hydroelectric power to lessen the cost of water deliveries, q 
water for irrigation, municipal and industrial use, recreation, and fish and © 
wildlife will all be included in this one great project. The major c costs al 
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FRESH WATER 


Valley, San Luis ‘and Santa Barbara counties, and Southern 
California from Los Angeles to San Diego. Even if fresh water from the sea 
_ were available for municipal use in southern California, its use would defeat 
the irrigation features of this project in the San Joaquin Valley, as irrigation 
would be uneconomical without municipal use to share in the aqueduct og 


.* 1990, when additional water will be required, a new economic analysis can 

be made of costs of fresh water fromthe sea compared to the long aqueducts, 
_ before a second state aqueduct is constructed. 
i Early Municipal Water Aqueducts. — This California Water Plan is the out- 
growth of such great earlier works as: -1905- 13, the Los Angeles-Owens © 
River Aqueduct, 240 miles long, extended in 1930-40 and additional 100 miles © 
to Mono Basin; 1913-34, the Hetch-Hetchy Project for San Francisco, 167 
. a miles; ; 1925-31, the "Mokelumne Project for the East | Bay Municipal Utility 
District, 100 miles; 1932-41, the Colorado River Aqueduct of The Metropolitan 
—— of Southern California, 300 miles; Denver’s 32-mile tunnel 


through the Continental Divide; and many other notable municipal water sup 
Major Irrigation Projects.—In the field of irrigation, flood control, and 
_ hydroelectric ‘power generation, great works of the West constructed by the 
United States are too numerous to mention more than a few major projects, 
such as: Hoover Dam and the All-American Canal; Grand Coulee Dam and the 7 


Columbia Basin; Shasta, Friant and Trinity Dams and long canals of the = a 
_ California Central Valley Project; the Colorado River-Big Thompson Trans- 
Ble Divide Project; the Glen Canyon-Upper Colorado River Storage : 
| Project ; and the great works of the Corps of Engineers and the Bureau of 
Reclamation onthe Missouri River, 


John Wesley Powell, who first navigated and | explored the Colorado River 4 
ont gave an intelligent analysis of the opportunity for man in the development © 
of the West. He had the vision that is still being realized and unfolded in the © : 


~ semiarid west. It was his ideas, enthusiasm, and work that led to the founding 
of the United States Geological Survey and the Reclamation Service, now the 
United States Bureau of Reclamation, Dept. of Interior, (USBR). _ The lake be- — 
hind | Glen Canyon Dam is to be named in his honor, as Davis Dam on the 
: Colorado River was named in honor of his nephew, Arthur Powell Davis, the - : 
early director of the Reclamation Service. 


‘These important events and outstanding works are mentioned to have 


the imaginative and zealous manner in which the settlers of the West have 
been undaunted in the magnificence of their works to capture water and bring 
it to lands and people. ‘The vision and deeds of the pioneers are succeeded by ~ 


enthusiasm of the West. In further support is the inherent faith of western 
people in the ability of modern science to find an answer, as it has so con- 
i spicuously | done in finding © the secrets of nature by research, An example of 
the results of the faith in, and support of, science in California is that nowhere _ 
: in America is there such an aggregation of Nobel laureates in science as at 7 
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™ located to water use are to be fully returned to the state by those benefited — ; 
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the California Institute of Technology, the University o of California, and ‘Stan- 


The low average rainfall and runoff of ‘the s¢ semiarid West in comparison to. 
that in the humid East is shown in Table 2. dl eo 


= TABLE 2.—AVERAGE RUNOFF AND WATER USE IN T 


Western | 3 


(millions) 
Runoff in Acre-ft 

(millions) 


of 
Population, 1950 
1960 
(millions) 


10 


‘drainage areas of the West, ru runoff off depths are as follows: 


Drainage Area ‘ Runoff Depth, ‘in n inches 


‘Missouri River - Bay Fae 
Rio Grande - Western Gulf 
Arkansas-' -White- Red Rivers 
‘California-Pacific 
‘Pacific Northwest 
_ Average of Western States 


These figures tend to indicate that “Pacific coast states es have 
q more favorable situation than | others of the West, but the remaining f fourteen — 
western states may expect severe limitations on their full development, -espec- 
jally for irrigation, But even in these three states there is severe maldistribu- 
tion of available water both in time and location. There is insufficient water 


full in the Snake River ‘Valley, eastern Oregon, in California 
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FRESH WATER 
generally so 
7 drainage areas. Exhaustive studies have demonstrated that there is sufficient 
= in California, Ancluding i its Colorado River diversions, to fully a 
ultimate requirements of some 50 million population and nearly 20,000,000 
irrigated acres. The necessary large storages, surface and underground, o 
- long aqueducts, hundreds of miles long, combined with high pump lifts, indi- : 
that the costs will be high. 
Need for more conserved water in the semiarid West ‘(with the 
_ favored locations in the Northwest and northern California) is a truism that is 
' recognized almost universally. Westerners have been conditioned to supporting | 
large expenditures of their savings and means to construct works to assure 
them of / adequate water supplies. They have become accustomed to look to the 
leading engineers to find the most economic solution to their future water 


~ needs. It is with this thought and background that this paper is written. 
MORE LAND THAN WATER IN THE — 7 
There is far more land than can ever be irrigated ne ‘the streams and 
‘ground water r of this, on the whole, _water-deficient area. There are i. 500° 
; miles of lands facing the sea. The Office of Saline Water shows that there are 
7 320 x million cubic miles of sea-water. For all practical purposes the quantity — 
is unlimited. There are also large amounts of inland saline waters. In fact, 7 
return water from irrigated lands is contributing to these all the time, and 


- their quantity is increasing. Almost all the — s saline waters, with the 


moteness from sea-water and difficulty in ‘disposal of wastes from inland 
of fresh water from saline waters. 7 ats 
_ What are the needs for developing fresh water from the sea along the Pacific 
coast of the United States? In Washington and Oregon, there is a great surplus — 
“4 of. excellent fresh water, , making the use of converted sea-water permanently | 
mye The same situation is true of the California coast north of Marin 


County in the ‘San Francisco Bay area, ', where there | are great surpluses of low- a 


where the comparable costs of water recovered from the sea and natural water 
transported in long aqueducts from the Sacramento- ‘San — Delta to 
water- areas must be, and have reviewed. 


CALIFORNIA VOTES 


CALIFORNIA VOTES $1.75 BILLION FOR INITIAL 
STATE WATER PROGRAM 


a are to be found in the West. There are Frond problems of elevation and re- 


Fortunately, 1947 has been engaged ina thorough study 
available water supplies, present and future water requirements, and the 

‘most economic means of conveying surplus waters from the north to meet the 
these: supplemental water needs of the balance of the state. The results of a 


these studies: are available. 3 They clearly ir indicate that there is sufficient 
natural water, including | Colorado River water, to meet all future water needs 
within the state. Large surface ; and eaenee water storage facilities and long 


a ueducts will ‘be re uired. 
= 3, State of California Water ‘Resources, , Calif. = lif, 
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The State Water Resources Development Bond Act to $1, 750° 
- ‘million in bonds was approved at the general election, November 8, 1960. 
_ ‘The average cost of water delivered wholesale to southern California is esti- 

mated to be $50 to $70 per acre-ft. It is planned that these works will meet 4 f 
the water needs” of definite service areas which contract with the state from 
+ the north San Francisco Bay area to San Diego. Costs per acre-foot will <a 7 

-much less for shorter distances and lesser pump lifts along the aqueduct. a 

WATER SEA CANNOT NOW COMPETE WITH 


While it is true ue that : for some small isolated cases, i. as that of the in- 
land oil town of Coalinga, a water supply from local saline water may be the 
present answer, the high production cost of fresh water from the sea for large 7 
f quantities cannot compete with the long aqueducts in California where iver, 
of gallons of water per day (millions of acre-ft per yr) will be delivered. 
‘The cost of the water at the Delta under the California state multiple-_ 
purpose program is “only 1.075 cents per 1,000 gal ($3.50 per acre-ft). For 
vi use in the San | Joaquin Valley, fresh water from the sea would have to compete - 
with this: cost. All other costs of the state’ program are the costs of f pumping» 


ors Going now to the "other areas of epportentty which may exist for fresh 
oa. extracted from sea-water along the Gulf Coast, Louisiana may be dis- 
“missed because of having access to large quantities - of fresh water except, 


: ‘becomes less to the west and ; south, diminishing to approximately 25 in. at 
Brownsville. coastal plain of Texas differs from that of the Pacific coast. 
Texas has slowly rising lands continuing 100 or more miles inland from 7 7 

a Gulf, while all along the Pacific coast the Coast Range of mountains parallels 
the coast, leaving but an occasional depth of plains or low, rolling hills for _ 
large population, industry, and irrigated agriculture. The largest of these 
ee plains is occupied by the Los Angeles metropolitan area. _ Pe 

 —s- The total | runoff of Texas rivers and streams is much less than that of 
California, and usable lands are several times larger. . Accordingly, it is con- 
ceivable that there may come a time when this region may look to the sea as 
an additional source of fresh water, especially for industrial use, This is, 
of course, , dependent on delivery « of such water at reasonable prices. a appears — 


the time is far off. 
PROBLEMS OF PRODUCING FRESH WATER 
INLAND SALINE WATERS 


- the saline waters of ‘the country. In general, such waters do not exist in the 
- Columbia River and coastal drainages of the Pacific coast. They do exist 
_ throughout t! the Rocky Mountain states , and are extensive _ through the Great _ 
Plains from the Dakotas to Texas. One ofthe problems of continued irrigation 
of lands is necessity of washing the salts downward so they will not 
“accumulate ine the soil. This: drainage frequently a serious 
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isposal of Concentrated Saline Water Wastes.- —A low-cost means of pro- 


benefit. . Its extensive use, however, is subject to” careful analysis of effect | 
on streams and ground water basins. In any process for the production of 
fresh water from saline waters there is the necessity of getting rid of the 
more concentrated saline water thus created by extraction of the fresh 
_— In all processes there is no reduction in the total volume of salt, ‘and : 
Be. every process having a reasonable prospect of ultimate use, it is not dry 
galt that remains” to be disposed of, but salt solutions of of varying concentra- 
"Sea and the Great Salt Lake andagreat many lesser surface and ground water 
accumulations of saline water, would have no great difficulty in disposing 
7 of the more concentrated saline water. Through the passing of time, however, 
the available raw saline water would become more and more en, 
adding to the costs of recovery of fresh water. 
2 Many states have already encountered the troublesome p problem of disposing 
of oil field brines, which may vary from 5% to 95% of the volume of oil ex- 
tracted from the oil- bearing formations. ‘The release of such brines, which — 
may” ‘be much more concentrated than sea-water, has severely polluted many 

: inland streams and ground water basins. Frequently, the only remedy is to 
oo the waste brines from oil production back into such deep porous rocks | 
as those from which the oil with water was initially extracted. This is by no 

_ ‘means inexpensive, though sometimes it has other offsetting values in in- 
creased production of oilorgas. 
An example4 of the high cost of pumping waste brines into deep under- 
_ ground rock formations where they will not contaminate useful fresh ground © 
SB canagl or surface streams is givenby these figures: 8,948,763 barrels of brine 
were so pumped underground at a stated cost of $348,926.94, or $0.039 per 
: _ barrel. _ This may not be a high cost in relation to oil at $2 to $3 per barrel, 7 
- ‘but converted to a water cost of $0.93 per 1,000 gal ($03 per acre-ft), it would 
be a heavy financial burden on any fresh- water plant which has to Gapees of 


its waste brine inthis manner. 
ie © CONCLUSIONS 
it is not enough to emphasize the increasing demands on our fresh-water | 
and suggest that future shortages, when they occur, may be met by 
salvage of fresh water from the sea and inland saline waters. ae See 
Insofar as the United States mainland is concerned, it would appear that | : 
there will be no substantial need or opportunity for such plants along the Paci- - 
tie coast. A. need for these plants n may develop along the Texas Gulf coast when 
the lower-cost natural waters are fully utilized. Inland in the severnteen 


western states, there will develop a great need for additional sources (of water, 


and there are substantial quantities of saline waters. This is an area where 
agriculture, especially irrigated agriculture, has a prominent or predominant 
— “Fole fa the present economy of each state. As natural fresh waters become 
cin developed it will be difficult for agriculture on a price basis to compete 


= “Oil Well Brine Problem of the Corney Drainage System (Arkansas- Louisiana), ” 
by Murray Stein, Twelfth Purdue Industrial Waste Conf., Purdue Univ., May 15, 1957. 
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4 with urban _ industrial demands for water. Any solution offered by production 
a - of fresh water from saline waters» would have to be within the ability of the 
= to pay. It would appear that high-cost disposal of waste saline- 
process waters in many instances, omer added problem in the use of 
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These problems are mentioned, because, int 

i should parallel that of developing lower cost methods of developing fresh 

water from the seas and inland saline 


in 


a A solution to the ata differential equation describing unsteady flow 
‘toward equally spaced drains above a boundary is 
presented, 

_ The solution is compared with field data ns a one numerical solution. 


On agricultural lands \ underlain of zero slope 
or where there is little natural subsurface drainage, equally spaced drains 
» may be installed to adequately lower the water table for crop production. 


The depth | and spacing of the drains are two important design factors that a 
“control the lowering of the water table. 


"ace open until Nov. 1, 1961. To extend the closing date one month, a 7 
written request must be filed with the Executive Secretary, | ASCE. This paper is part | 
of the copyrighted Journal of the Irrigationand Drainage Division, aati al « the 


; ~ American Society of Civil Engineers, Vol. 87, No. IR 2, June, 1961. 
Agric. Engr., Ft. Collins, Colo. 
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7 The solution was found to agree with the field data and the numerical solution Ee 


Proceedings of the American Society of Civil Engineers 
4 TILE 
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‘cent years, some senieamacail 3,4,5,6 have been obtained for the 


: ium of drainage problems involving a falling | water table. Many of the 


approximations are objectionable? because the assumptions on which they are 
oe, based are not realistic and hence, the solutions may not be sufficiently accur- - 


ate for design purposes. . This discussion deals with improving an approxima-_ 


tion developed by R. E. Glover, 3 F . ASCE, by eliminating an assumption . _ 


his development. Such an approximation will provide information < on — and 


As reported by L. E. Dumm3 Glover developeda solution on the 
flow equation for the problem of the falling water table. Using the Dupuit © 
assumption, he considered equally spaced tile drains in a homogeneous soil 7 


overlying a horizontal impermeable layer. In development of Glover’ 


solution an approximation was made that restricts the use of his equation to 
= cases” in which the distance from the tile to the impermeable boundary is 
large in relation ‘to the | drawdown of the water table. Experience has indi- 
that when the “assumptions are Satisfied, satisfactory results may 
obtained. In summarizing some-of the recent studies on the falling water table, 
horizontal flow, appears | to be. ‘more nearly correct | than any other, but it is 
not sufficiently accurate to be used for design purposes.” = 
me _ The approximate solution presented herein takes into account the drawdown 
of the water table which implies that the distance from the tile to the imperme- 
able ‘boundary may be small compared to the drawdown of the water table. 


van Schilfgaarde? ‘stated that “Glover’s equation, based on the assumption 


for flow ba: based o on the mn the Dupuit assumption and Darcy’s s law may be written as ~ 


Model Tests on a Tile-spacing Formula,” by William Ww. Donnan, Proceedings, 


Soil Science Soc. of Amer. , Vol. 11, 1946, p. 


3s “Drain Spacing Formula,” by L. E. Dumm, Agricultural Engineer, Vol. 35, 1954, p. 


4 “A Quantitative Method for Determining Ground Water Characteristics for Drain | 

, Agricultural Engineering, Vol. 31, 1950, p. 284. 

- Landwirtschalfliche bodenverbesserun, by J. Spottle, gen Hamb. Ing. Wiss., Part 
_ 3, Wasserbau, Vol. 7, 4th edition, Wilhelm, Englemann, Leipzig, p. eo, ou vail 
«6 “Depth and Spacing for Drain Laterals as Computed from Core Sample Permeabil- 


“ity Measurements,” by P. Walker, ‘ Agricultural Engineering, Vol. 33, 1952, p. 71. ee 
_ 7 “Approximate Solutions to Drainage Flow Problems,” by Jan van Schilfgaarde, 
"Drainage of Agricultural Lands, edited by James N. Luthin, Agronomy, Vol. 7, 1957, 
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‘lis ‘substituted into of ‘continuity, 
flow toward the 


the « diffe rential 
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in which D f is the yield, and t denotes time. 
he and initial condit on the system are given as 
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Emile” ‘Picard’ 's method of successive approximations® wa was s suggested to 
___ SOIL_SURFACE 


linear differential Eq. ‘2. Eq. 2 is not 


{IMPERMEABLE _BOUNDARY 
2 


-MODE L OF LLY SPACED DRAINS ABOVE AN IMPERMEABLE BOUND- D- 


exactly the 1 same type of equation considered = Picard, however, Picard’s - 


method when applied to this equation yields the same results as the method 
‘of Poincare, Lighthill and Kuo, sometimes referred to. as the PLK method. r. 
The application of the PLK method to Eq. 2 has been clearly presented by 


9 “The Poincare- Lighthill- Kuo Method,” H. ‘Tsien, Advances in Applied Me- 
chanics, Vol. IV, Academic Press, 1956, p. 281. 281. 
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_ 8 “Mémoire sur la theorie des equations aux derives patielles et la methode des 
_ approximations successives,” by Emile Picard, Liouville’s Journal de Mathematiques, 
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x 


in € is the perturbation parameter. When € = 0, reduces 
linear differential equation that was solved by Glover as a first soproximation 


= Assuming that the solution h(x, ‘ € ) may be expanded in in ‘powers « of €, ‘the — 


h(x, t, €) 


substituting Eq. 5 i into 4 and equating the coefficients of of € 
to zero, a set of differential equations is obtained, the solutions of which © 


F provide the necessary terms in Eq. 5. This set of differential equations is 


‘a 


4 
The first teami in Eq. 5, ho(x, t), is Glover’s solution and oo it is modified 
for the boundary conditions given by Eq. 3 becomes peat 


_an cor 


10 Discussion by K. T. Raja Iyengar Sundara of “Unsteady Flow of Ground Water into 
a Surface Reservoir,” by William Haushild and Gordon Kruse, Proceedings, ‘ASCE, Vol. 
87, No. HY 1, January, 1961,p.280. 
“Unsteady Flow of Ground Water into a Surface by William Haush Haushild 


i, July, 1960, Pe aed 


— 
— wherein Eq. 2 was dierent Set of boundary 
As stated by S 
— &g 
— 
| 
= 
= _ The second term in Eq. 5 was obtained from the solution of the differential _ : 
- equation in Eq. 6 involving hy andho. Each solution for the Eqs. 6 must satisfy = J 
| 


the conditions given by Eq of Because of the difficulty in obtaining additional 
terms, the solution presented subsequently involves ; only the first two terms 


of Eq. 5. The solution to the problem is — in terms of hy &, t) and is 


m= m= L 


= 
term, G, under the time is equal to to + 52) that 


. function of [(t - T), x] in : which Tis a time variable that takes on ote s range ve 
values t > T> > 0. The integral 
- not satisfy the ne vation, hence, the time integral in Eq. 8, that 
was evaluated numerically, introduced to permit restoration of t the 
A solution similar to that p presented by Haushild and Kruse is also | pre- 
g sented here for purposes of comparing solutions. A peeeecuael approximation 
for flow toward the drain in Fig. 1 is given by t- 


* in , h < <) ‘D. A more accurate expression for flow toward the drain is. 
given by Eq. 1 


oe 


An approximate solution, hj, is | 
when integrated yields 


ion, Using th the boundary conditions given 


by Ea, 43, the constant of integration is found to be/ -After substituting the 


7 
A(x, t) = (1 
(| 
| 
i 
» 
11 
y in which c is the constant of integrat aie 7 
| 


e of c into Eq. 12 and rearranging, the new vequecstmetion, hi, is expressed 


cua The solution of Eq. 8 is presented graphically in Fig. 2 for x = 0 and for 
“various values of H,/D. Glover’s solutions are also shown in Fig. 2 for com- . 
parison with this theory, 
Behe The reason for the apparent inconsistency between Glover’s special equation 
for H,/D = 2.0, that places the drain on the impermeable boundary, and Eq. 8 
is that Glover’ ’s solution? does not satisfy the same initial conditions assumed ~ 
i ‘The « curves in Fig. 2 obtained from Eq. 8, ‘show that the water table midway 
between drains does not begin to recede untila finite period of time has elapsed 
_ depending on the distance from the drain to the impermeable boundary. The - 
range : of relative distances to the impermeable boundary, 0 <H)/D < 2.0, con- 
stitutes 100% of the total drainable depth. Don Kirkham!2 showed that flow into | 
tile drains is considerably reduced when the drain is placed on the imperme- 
_ able boundary. A comparison of the curves in Fig. 2 for the drain on the im- - 
‘ permeable boundary, H)/D = 2.0, and at two-thirds of the drainable depth, 
_ Ho/D = 1.0, also indicates that when the drain is near the boundary the flow is 
comparison of the two approximations, Eq. 8 and 13, is shown in Table 1. 
oo agreement of the two approximations is excellent for 0.1< H,/D< 1 1. a 
_ However, for H)/D = 2.0, the agreement is not entirely satisfactory, 


ee a recen paper by erwoo e general method o n Kirkham 
he: In t by J.D. Ish d13, th 1 method of Don Kirkh : 


and a Gaskelll4 was used with the use of a high speed computer to solve — 
the boundary value problem presented herein. Isherwood solved the problem 
as a two-dimensional system; hence, his solution should be more accurate than = 
_the one-dimensional system which is based on the Dupuit assumption. By com-~ 
paring Isherwood’s solution with the approximate theoretical solution pre- 


sented here, it is possible to de determine the validity of ‘Eq. 8 to a reasonable ~ 
— 


; A Isherwood’ s data for the case of an impermeable boundary existing 10 ft 


& 


below the ‘soil surface is reproduced in Fig. 3, , using dimensionless para- 
a restriction 


in in Soil Underdrains Resulting from Their Partial 
"ment in, or Proximity to, an Impervious Substratum,” by Don Kirkham, Proves’? a 
Soil Science Soc. of Amer., Vol. 12, 1948, p. 54.00 
a “Water Table Recession in Tile-drained Land,” by J. D. Isherwood, Journal of 

Geophysical Research, Vol. 64, 1959, p. 795. a 


so "Falling Water Table in Tile and Ditch Drainage,” by D. Kirkham and R. E. -_ 


kell, Soil Science Soc. of 15, 1951, 37. 
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(and Eq. and This restriction may be made by use Of the ratio 

1 9 ile spacing, L, to the distance from the tile to the impermeable boundary, | 
P 
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FIG. 2. 2.—- THEORETICAL , CURVES OF RELATIVE WATER TABLE HEIGHT, h/Ho, 

A FUNCTION OF THE TIME PARAMETER qatt/L2 FOR TILE AT VARIOUS | 

RELATIVE DISTANCES Ho/D ABOVE THE IMPERMEABLE BOUNDARY AND 
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i “that is, L/d. It is reasonable to expect that for large ratios of of L/d a solu- 
mate ‘the ‘two. dimensional solution, and conversely for ‘small ratios of L/d 
such a solution might be in considerable error. This fact is clearly shown in 
 «*Fig . 3 in which Isherwood’s solution agrees closely with the theoretical curve 
for L/d >25 whereas for L/d< 25 there is considerable lack of agreement. It p 
should be noted in Fig. 3(a) that Isherwood’s data for L/d = 25 coincides with — 
the theoretical curve of this paper. In Fig. 3(b), Isherwood’s data for L/d = 50 
— is represented as points because the two curves for L/d = 25 and L/d = 50 are 
; nearly coincident. It appears from the data available from Isherwood’s paper _ 4 
: a: that th the e theory of this this paper is is independent of L/d for ratios greater than 25 


endent of L/d 


-TABLE 1.—A COMPARISON OF THEORETICAL SOLUTIONS OBTAINED FROM EQS. _ 
= AND 13 FOR = 0 AND FOR VARIOUS VALUES OF H H o/D. 


0.67 1, 


9 0,50 


3 
8 

15 

01 

.23 

31 
8 
3 


-0.4 


of F. -Klinge’s unpublished field data, and ‘Eq. 8 is 
perms in Fig. 4. The theoretical curve appears as a solid line whereas the 


theoretical cuve. “These « data were reported for shallow drainage systems in 
which the impermeable boundary was 3 ft to 4 ft below the surface. 5 — . 
of these comparisons were made and similar agreement was found. = “a 
a An example of the use of the theoretical curves follows. ‘Assume a system 
: ot drains are to be installed to a depth of 6 ft below the soil surface and ~_ 
impermeable boundary exists at 9 ft below the soil surface. If the position of 
the water table initially exists at the soil surface, compute the tile spacing 
_ that will cause the water table to recede to a depth of 3.6 ftin3 days when the — 
hydraulic conductivity-specific yield ratio, K/f, is 2.0 ft per hr. The scaled 
time variable, t/L2, is read from the theoretical curve, H,/D = 1.0 ) for 
h/ ‘Hy =- 9. '-0.1, that is found t to be 0.150, Fig. 2. Because t and a are known, , the 


a, 15 *Physical Giameinamtdalen of Clermont Silt Loam Soil in Relation to Tile Drain- 


age,” by A. F. Klinge. Thesis presented to Purdue Univ. at Lafayette, Indiana, in 1955, = 
in partial fulfillment of the for the of Science. 
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in which a is 12.0 ft2 per hr. The spacing is found to be 76 ft (approximately) © 
and L/d = 76/3 = 25.3 that is within the limits set forth using the data : 

Isherwood. The theoretical predictions of the tile spacings that can be 

tained from Fig. 2 willbe nobetter than the methods used to obtain character- 
istic hydraulic conductivity - specific yield data for the drainage system. 
_ Methods of measuring hydraulic conductivity on a large scale are needed. - 

_recent: work Of seems to be a step forward in this di- 


“Will demand the study of the soil as 
_ both above and below / the water table,as an essay in the field of water econ 


— 


- ment in a medium whose hydraulic conductivity is a function of moisture con- 
é tent. ” However, it is hoped that the approximate solution presented herein will 


_ enable the drainage engineer to more accurately design drainage installations — 
predict drainage costs. 
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WATER TA FLUCTUATIONS INDUC ED BY 


Marinus Maasland, 1) M. “ASCE 


“as a he theory of intermittent recharge, as previously developed, is extended 
and analyzed. Equations are presented from which charts may be prepared 
-_ _ for easy computation of drain spacings for irrigated agriculture. The analytical 
results| may also- be used to compute the effect of natural drainage onthe 
- height of the water table as well as anticipated changes in the water table of _ 
newly” irrigated areas. Example applications of the theory to ditch and all 
and natural drainage will be covered in later publications. 


The effect of recharge on the table on many tactors 
characteristics of the water bearing stratum and of the materials ‘under and — 
% overlying it; amount and distribution of the recharge; extent to which the river | 
or drain penetrates into the water transmitting porous medium; and the depth = i: 
spacing of subsurface drains. The theory of intermittent recharge2 has 
been found useful in analyzing relationships between the various factors. The _ 
mathematical treatment is somewhat complex but its results can be — 


iently used for solving practical problems by preparing graphs ; and charts. 


analysis brings together, in a rational manner, , the opposite 


- Note. —Discussion open until November 1, 1961. Toextend the closing date one month, 
; a ‘written request must be filed with the Executive Secretary, ASCE. This paper is part 7 
e of the copyrighted Journal of the Irrigation and Drainage Division, Proceedings of the - 
Society of Civil Engineers, Vol. 87, No. IR 2, June, 1961. — > 
; 1 Engr., Tippetts- -Abbett-McCarthy-Stratton, Engrs. and Archts., 375 Park Ave. —_ 
_ 2 “Watertable Fluctuations Induced by Intermittent Recharge,” by M. Maasland, Jour- 
of Geophysical Vol. 64, ~559; (a) Pp. 551; (b) 557; (c) p.5 556. 
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assumptions of instantaneous recharge? ar of 7 constent re 


_ charge. 2.45 Th The analytical results 8 may ‘be used | to determine the following: 


1. Drain spacings for irrigated @ 
= 2. Applicability and limitations of the Glover ‘formula’ an and of steady state” 


-4 equations for computation of drain spacings; = ~ ier er 
3. Combined effect of excess irrigation papumienng rainfall, and natural 


4. "Water table levels that ultimately. will be reached in aquifers under the 
7 5. Number of years needed 1 for the water table to approach the new equi- — 
librium level in newly developed irrigation projects, following increased ae 
Pe 6, Effect of various amounts and methods of application of excess irri- | 
water on the height and fluctuation of the water table; 
* Type and distribution of recharge received by the groundwater res reser- 
voir from detailed field studies of water tables wepeenienets of natural 


percolation losses irrigation projects. 


The equations aid in recharge conditions of ground- 
= water basins and drainage needs for irrigated areas. They also permit esti- ~’ 
"7 mating the time when drainage construction will be necessary for newly — a 
te veloped irrigated areas so that funds and equipment for this purpose canbe _ = 
scheduled. High water tables and attendant salinity problems may thus be 
by acarefully planned program. = 
; Tt is known that the first Glover formula3 ~ may lead to narrow spacings 
that are inconsistent with field experience. In applying that formula, the mis- - 
leading assumption is made that itis necessary to remove from the root zone, | 
or from the soil layers below it, the total amount of deep percolation loss from 7 
each irrigation prior to applying the next irrigation. Maasland indicates? that 


- following a single lueetinnenh recharge. It does not account for either the 
~ accumulative | effect of successive recharges or the depletion of the ground- 
water between or following a number of recharges. However, 
_ storage during the irrigation season and its subsequent depletion during ~ 
_ non-irrigation season have a significant effect on the water table level. By ~ 
: = these factors into account, it can be shown that permissible drain 
spacings may become much wider than those computed from the Glover form- 
ula. It will be shown elsewhere that the various methods of determining drain 
a and the spacing formulas derived, recommended, or used by William a 


w. ‘Donnan,® F . ASCE, et al., . D. Dumm, ‘Isherwood,’ 


ai 3 “Drain Spacing Formula,” by L. Agricultural Vol. 35, 1954, 
Bijdragen: tot de kennis van enige ge natuurkundige grootheden van de | grond, "by Ss. B. 
Verslag, Landbouwk, Onderzoek, Vol. 46, 1940, pp. 515-701. 
5 “Seepage of Steady Rainfall Through Soil into Drains,” by D. Kirkham, Trans- 
actions, Amer. Geophysical Union, Vol. 39, 1958, pp. 892-908 (a) p. 900, Eqs. 84- a7. 
 ——_ "Drainage Investigations in Imperial Valley, California, 1941-1951 (A 10-yr sum- 
_ mary),” by W. W. Donnan, G. B. Bradshaw, and H. F. Blaney, U.S. Dept. of Agric., Soil 
Conservation Service, Tech. Pub. 120, 1954. 7 
7 “Watertable Recession in Tile-Drained Land,” by J. D. Isherwood, Journal of Geo- 
physical Research, ‘Vol. 64, 1959, pp. 195- 804. 


| 

i 

g 
— 

= 

> 

| 

| 
i 
= 

® 

— 


WATER TABLE 
houdt,4 4 Maasland, 8 Maasland and H. C. Haskew,9J.D. Van Schilfgaarde,10,11 
et al., and D. "Kirkham often require adjustment of the recharge coefficient — 4 - 
4 they are to be used to compute drain spacings for irrigated agriculture ; 
a= errors may result if the recharge conditions are not properly consid- 
- The equations for intermittent recharge are derived on the cen 
that the drains completely penetrate the permeable stratum, If the drains a 
_ penetrate only part of the permeable stratum, Eqs. 29, 70, 77, 78, 82, and 83 
Ss Hooghoudt4 may be used to account for flow convergence near the drain, 8a 8a a 
The work of Kirkham54 may be similarly used and leads to approximately the 
~ results as do Hooghoudt’s equations. As a matter of interest, it is ob- 
a served that Hooghoudt’s Eq. 29 and Kirkham’s Eq. 87 are identical. This” 
identity follows immediately from Wallis’ infinite product for 77/2. Correction 
for flow convergence may be done by a graphical | method as previously pre- : 
® it is emphasized that this paper is restricted to an examination and ex- 
tension of the theory of intermittent recharge. Applications will follow in 
4 subsequent works in which appropriate graphs and charts will be presented, 
The charts, , which relate various dimensionless parameters, allow rapid — 


eva aluation of a of conditions. = 
r are defined a 


alphabetically, for convenience of reference, in the 


tig General. and ‘elaboration of the theory follows to to ex- 
_ ploration of the practical implications of the analytical results. Frequent ref-_ 
erence will be necessary to the writer’s previous work2 and to equations in 
‘it. For reasons of convenience and to prevent repetition, the previous paper?» 


will be referred to as (MM). Its equations will be indicated by their number — 


aa M. 33. The notation used ‘in the present paper is generally the same as that 
used in (MM). A few small changes have been made which will be indicated. 
> the previous paper,2 two related types of intermittent recharge were 
considered Separately; intermittent constant recharge (recharge applied - 
at a ‘constant rate), and intermittent instantaneous recharge 
(recharge applied instantaneously at regular intervals). The effect of inter-— 
_ mittent recharge on the water table and on the rate of discharge of the drain 
= (base flow’ in the case of natural streams, or ‘return flow’; a term used to 
denote the discharge induced by deep percolation losses from excess irrigation 


- applications) may each be studied by considering simultaneously two functions _— 


8 “The Relationship Between Permeability and the Discharge, Depth, and Spacing of : 
® Drains,” by M. Maasland, Groundwater and Drainage Series, Water Constr. = 
ee Pty Committee, N. S. W., Australia, Bulletin 1, 1956, p. 35; (A) pp. 21-33. at lle 
es & 9 “The Auger Hole Method of Measuring the Hydraulic Conductivity of Soil and its 
_ Application to Tile Drainage Problems,” by M. Maasland, Third Congress, Internatl. . fi 
4 Committee on Irrig. and Drainage, Question 8:8.69-8.114, 1957. = = © > We 
10 “Physical and Mathematical Theories of Tile and Ditch Drainage and Their Use- 
"fulness in Design,” by J. Van Schilfgaarde, D. Kirkham, and R. K. Frevert, Iowa State 
College Agric. Experiment Sta., Bulletin 436, 1956, pp. 667-704. 
4 11 “Drainage of Agricultural Lands, Agron,” by J.N.  Luthin, Amer. Soc. Agron., 
Monograph N No. 7, Wis., 1957, 620. 
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such as a - periodic one and a transient one. The height of the water table or 
_ the rate of discharge at any time is found by subtracting the transient from 
the: periodic function. The transient is identical to the periodic function = : 
ae _~ = 0, tbeing the time since the beginning of the first recharge. After sufficient: : 
time has elapsed the transient becomes negligibly small. The water table and 
_ the discharge are then completely defined by the periodic function. oo. 
The factor uw of (MM) was not given ‘a subscript n because it occurred so 
3 frequently in the analysis. The subscript was not needed in the derivations 28 
_ long as the dependence of uw on the number n of the term in the series was _ &§ 
_ realized. The subscript will now be used consistently because it is cco 


3 numerical value of the dimensionless 


of UoT (less than 0. the of fluctuation as caused 

recharge (of rate p= a/t) is small in ‘comparison to the. average water table 
- level or discharge. In that case, the effect of intermittent recharge is, ae a 
practical purposes, _ equivalent to that of : a continuous constant recharge at a — 

_ lower rate p; = q/T even if the interval of recharge et; is considerably less 
than T. The ‘simple steady state Eq. M.3 and the nonstationary Eq. M.7 for © 

continuous constant recharge may then be applied to problems actually 


volving intermittent recharge. The; amplitude of fluctuation can be estimated 
from q/e or a fraction thereof. It may be shown that, in the limit for UoT — 0, 4 
_ Eq. M.23 with p = q/t; is identical to Eq. M.3 with p; = q/T. The phenomenon ~ 
referred to is independent of the interval of recharge ty. Proof will be given 


later. For — values of Seales the pectiom of intermittent recharge degenerates 


in which hx 1 represents the height o of the water table measured from the line _ 

connecting the water levels in the drains. Note that hy, as defined here, equals — 

2 Ss - hy) of (MM). Eq. 2 implies the assumption that Eq. M.2 is applicable to 
all problems considered. Since solutions of Eq. M.2are additive, hy may also _ 
be taken to be an added height to the water table prevailing before the be- 


ginning of the newly i introduced intermittent recharge. lace . 


a If H (Fig. 1) is small 2 Eqs. M.2 and M.3 may be inapplicable, so that Eqs. 
4 M.5 and M.6 must be used instead. This requires that all equations of (MM) | 
_ be transformed as indicated in (MM). Za In addition, the analytical procedures» 
7 _ developed for intermittent recharge may then not lead to sufficiently accurate 
_ results if D increases significantly as the water table rises. _ Variation of D 
causes to change, which invalidates Eqs. M.17, M. 18 and. M.26, , because 
— Eqs. M.14 and M.15 and S of Eq. M.25 no no longer 1 represent geometric series. 2a a 
must then be taken to a procedure involving the use 
of the (properly transformed) Eqs. M.13, M.14, M.15, and M.25. Atabular 
_ procedure for dealing with these problems has been n developed by W. N. Tapp 
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a: 
| 
ij 
7 
| 
ij 
| 


IR2 WATER TABLE 
and Dumm. Another tabular procedure is used by Glover, F. ASCE, and 

M. W. Bittinger.124 To keep this paper within bounds, the methods referred ‘7 
to will not be examined herein. The principle of superposition, /used in both 
the tabular and analytical methods, is not a new one. It has been applied to _ 
various ground water problems by J. Boussinesq,/3 J. H. Edelman,14 Jacob, 
Hooghoudt,4 Kirkham,!1 p, A. Kraijenhoff,15 P, W. Werner,17 and others, and 
is widely used in operational mathematics. 18,19 
dae The symbols hp and hpj, and h¢ and hej were used in (MM) to define the 
water table level for any x. These symbols will now be used exclusively in — 


‘reference to the level midway between the drains. The eqrenaneys symbols © 
= me water table level for any x are given the subscript x, that is, ip,X and 


hpi x, and ht.x and hej _x. _ Thus the symbols | are given a subscript for the 
eral problem, while,’ for the special case in which x = 1/2L, the symbol h 
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FIG. 1.—GI OUNDWATER INTO OVER A HORIZONTAL IMPERVIOUS | 


_ et “Transient Ground Water Hydraulics,” by R. E. Glover and M.W. -Bittinger, Colo- 
-rado State Univ., Fort Collins, Colo., No. CER59REGI6, 1956, p. 57; (a) pp. 43-44. > 
rs. 13 “Note sur le Mouvement des Eaux Souteraines. Essai sur la Theorie des Eaux > 
by J. Memoires, l’Academie des Sciences, 23, ‘France. 


Groundwater), thesis presented to Delft in 1947, in fulfillment of the require- 
ments for the degree of Doctor of Philosophy, esi‘ 
1 “A Study of Non-Steady Groundwater Flow a Reference to a Reservoir- a 
- Coefficient, ” by D. A. Van De Leur Kraijenhoff, De Ingenieur, Utrecht, Vol. 70, 1958, 
7 PPS “On Non-Artesian Groundwater Flow,” by P. W. Werner, Geofis, Pura Appl., Vol. | 
17 “Some Problems in Non- Artesian Groundwater Flow,” by P. W. Werner, ‘Trans-— 
actions, . Amer. Geophysical Union, Vol. , 38, 1957, pp. 511-518. 
“Modern Operational Mathematics in Engineering, by R. Churchill, ‘McGraw- 
Book Co., New York, 1948, p. 306.0 
“Conduction of Heat in Solids,” by H. s. Carslaw and J. 
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remains s without a subscript identifying the horizontal distance from the drain. 
‘There is no need to further define the rate of discharge of the drain Q, since 
the (identical) discharges at x = 0 or L are the only ones considered ome a 


= —_ 


ts 


¥, are variables and their ur unit is. equiv- 

a in time to the interval T T; y varies from 0 to 1, while é may assume any 

er . periodic function for intermittent constant recharge consists of two 
equations. _ These equations, which were combined into a single expression — 

7 (Eq. M. 23) to facilitate the analysis, will no now be stated separately f for the 
two intervals and their limits. | The equations» for the rate of discharge, 4 


_ instantaneous recharge will be examined in some detail, 


, derivation of which was indicated (Eq. M.39), are also stated. Intermittent 
| 
‘In order to evaluate the effect of drainage s systems on the water table for 


irrigation agriculture, the theory will be extended to include the problem of a 
‘double periodicity.’ | It is convenient to divide each year i into two such 


requires only a slight extension of the theory developed for ‘single 


‘ a restrictive one and makes the theory inapplicable for certain areas. It 7 
_ Can be applied only if the recharge during the nonirrigation season is ne 


- as compared to the recharge during the irrigation season. : - is, of course 

possible to extend the theory to include recharge during the nonirrigation 
i ‘However this would further complicate the equations and the writer 

plans to. show, in subsequent publications, that such complications are not 
usually necessary. If a highly irregular recharge distribution must be Anal 

‘sidered, it is desirable to use a tabular rather than an analytic procedure. A 

procedure for ‘tile: drain spacings has been described by 
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it will be assumed that the ground water is recharged at regular 


WATER TABLE 


bon The Water Table for Intermittent Constant Recharge. —Eq. M. 23 isthe gen. gen- 

eral expression for the periodic part of the equation for the water table. This | 
_ expression must be evaluated for particular values of 7 (in whichO S 7 s T) 
is applicable to any interval (r - - If the variables 5, 


‘defined by ‘Eqs. 3, are entered ir in 1 Eq. M.2 23, it becomes 


) exp - Hy Ty - exp - Hy To 
n which (expy. Tt -1 


Eq. 5 may be w written 

A 


7 - of Wx is is used to indicate that this expression 1 is valid for0<x < 
6 


th the case of intermittent constant recharge, each period T consists of | 
two intervals such as one during which recharge takes place (at a rate p 
= q/t1) and the other without recharge. Due to Eq. 4, the periodic function | 
4 consists of a separate equation for each of these intervals. These equations ¥ 
do not have the same form but are identical at the limits of the intervals, that 7 
is, for y = 0 and 1, and for y = ¢. Tt has been 1 emphasized (MM) that Eq. 5 is 7 
_single-valued for all values of y. 
a 4 For the special case that x = 1/2L, the equations for the intervals of y and 


«af 


Eq. 8(a) follows immediately from Eq. -9(a) by putting y = 0, or from Eq. 
10(a) by putting y= 1(MM). ForO<y<¢, Eqs.4and5yield 


= 


— 

te 

| 
or 
tm 
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= 


= 
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‘The maximum water table level, ‘defined ash occurs aty=¢ and is given 


by Eq. 9(a) or Eq. M.37. For ¢< y <1, it follows that : chet 
T 6 - An) exp | T (10a) 


, respectively, because 


‘the latter equations more than Eq. 10(a) for the ‘Par- 
x # the factors w of Eqs. 8(b), -9(b), and 10(b) may be wri as 
Why Wo x, and W3, => respectively, to conform with the notation previously » 
“bet for hp, x and ht. x x. The water table level for x # 1/2L is obtained © 


by substituting sin (2n + 1) t (x/L) for (-1)" in these equations. = a 


(3). sin (2n+4) (=)a, exp- Te 


which Té =t, the time since the beginning | of the first ‘recharge; r 
4 


The for the transient (Eq. M.24) is rewritten as 


assume any value from zero to infinite 


_ The Rate of Discharge for Intermittent Constant Recharge.—The rate _ 


per unit of drain (Q) may be from » using 
@ — co (tas. Ma) | 


hy 


= 
“The discharge Q is separated into its periodic and transient components _ 


Eq. +5 and the first of Eq. following expression is ‘obtained 


Although hp1 and hph fr 
ij 
- 
— 
7 
- in which Q) denotes the rate of discharge computed from the periodic function, a 
—— and Q, is the rate of discharge computed from the transient function. It is — 
. 7 noted that Q as defined herein represents the inflowfrom one side of the drain 
— in which d is that part of Eq. 14(a) between braces; d is determined entirely , 
by the dimensionless parameters T, y, andt. 


The minimum rate of discharge, defined as min» occurs at y = Oorl, 


and, _ from Eqs. 4 and 14(a), is given 


g\ 
(ones +1)" . (16a) 


2 T 


‘The rate of discharge as computed from the transient (Q) is obtained by 


iets the second part of Eq. 12 for == 0, using Eq. 11. The result is © 


Bari” exp-u, TE 
Tare 


én Eq. 18 is monotonically decreasing and becomes negligibly small after suf- : _ 
_ ficient time has elapsed. When this state of dynamic equilibrium is reached, 
<a the rate of discharge is determined entirely by the periodic function, and the — 
total discharge (Dr) per period T = will be equal t to the total recharge per 
- The variation in the rate of discharge decreases as Ho T decreases, and 
_ it may be shown, as is done for the water | table in a later section, that steady > 


state flow would obtain in the limit for Lo T— 0, after dynamic equilibrium — 


: has been r reached. It also follows that, for low values of Uy T, many periods 

: T are needed to approach the rate of discharge as computed from the peri- _ 

odic function, in other words, to reach dynamic equilibrium. = ~ 
1 


one, Se 
__ Intermittent . Instantaneous Recharge.—The equations for the water tab 


level for instantaneous recharge need not be repeated herein, as they have 
been adequately covered in the previous paper (MM). The water table level 
at any time is given by Eq. M. 28, _ while Eq. M.29 is the periodic function and _ 
Eq. M.30 is the transient. The minimum water table level for any x andfor ~ 
"dynamic equilibrium (Qpim,» x) is given by Eq. M.31, and the minimum equi- 
librium le level between the the drains (hpit) by Eq. M 32. levels 


are re obtained bya adding Eq. M.31 or Eq. M. 32. The symbols hpi, x and hej, 

as used herein refer to the water table level for any 3 x, , and d hpi a and hij are 
used to define the level between the drains. 


: 
| 
 For0sys 
| 
maximum discharge, de irs at y = ¢ and is computed 
(2 n4 - an) exp Lt, T ...(17a) 
Lore (1m) 
| 
| 
| 
| 
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ia The rate of discharge | per unit length of drain (Qi) may bec obtained from a 

The discharge > Qi is is separated into its oe: and transient components 


in which Qpi is the rate of discharge computed f: from the periodic function, and 


Qh denotes the rate of discharge computed from the transient function. Itis: a 
= that Q; represents the inflow from one side of the drain only. Ss ¢ 


‘From ‘rom Eq. M.29 and the first part of Eq. 19 it follows that 


. M.30 and the second part of Eq. 19 


— 


(23) 
‘It is remembered and that t t= T+ 7, or = (r-1) 
thal ‘Eq. 21 does not converge for 7 = y Physically, this is to be expected be- —_ 
cause instantaneous recharge ‘implies. the of recharge at 
rate. The rate of discharge is finite for 0 S 7 S T and maybe computed from — 
Eq. aged However, the Fate of convergence of the series in n Eqs. (21 and 22 is 


It is often more 
IP “a convenient to compute the volume of the water stored below the water table — 
7 and to determine the rate of discharge from the decline in the storage. The 
needed “computations may be shortened by computing for various times the — 
amount of ground water stored below the water table and plotting the ‘mass- 
-curve.2 7 The rate of discharge at any time is then obtained by drawing a 
- tangent to ) the mass- curve. This method may also be useful for intermittent 
_ constant recharge if Eqs. 14(a) and 18 converge slowly. A description war 
application of the principle of this method may be found in the work of Glover. ia 
and others. ate The procedure derives its merit from the fact that ee 
_ of the series in Eqs. M.29 and M.30 improves their rate of convergence, while 
_ differentiation has the opposite effect. Kraijenhoff considers, 15 in detail, the a 
storage may be determined as follows. ‘The water level for any 
x is given by Eq. _M. 28, and the volume of water stored below the water table — 
at any time is obtained by integrating Eq. M.28 with respect to x from 0 to L- 


cop 20 “Water Supply Engineering, ” by H. E. Babbitt and J. J. Doland, McGraw-Hill “oa 


21 “Cooling of Concrete Dams,” by R. E. Glover, Boulder Canyon Project, Final Re-_ 
ports, U. S. Bur. of Reclamation, Bulletin 3, p. 236; (a) 


— 
— 

a and from Eq g 

; 

> 

ine 

* 


multiplying the | result by €. 


in which S; denotes the volume of ground water stored below the winaiiie’ 
«per unit width of aquifer and between two: drains spaced at a distance L. ie le 
_- The volume §; will be separated into Spi and Sj, the storage computed from — 
the: periodic function and the transient, respectively. It is found from the first 
term of Eq. 24that : 
and from the second te term 


Te. (26) 
the fraction of the original volume of drainable 
water a+ in the aquifer. ‘The portion which has been discharged | through 
the drain is equal to( ay The mean water table is given by the expression — 


: 26 is monotonically decreasing negligibly small after suf- 


total discharge (Dr) per period T will be equal to the recharge per period, 
4 that is , 1/2q L. The factor 1/2 i is used herein aenageete the inflow from only 


-fetent time has elapsed. After dynamic equilibrium has been reached, the = 


to 7 which 7 = T) from 0 to or 


Using the 


itt follows immediately that 
14(a) for intermittent 
a _ The Factor a,. —The factor a, (Eq. 5) occurs in all problems of intermit- 
~ tent constant recharge. ‘The expression for a, maybe simplified for both high - 
and low values of Ho The simplified expressions may be used if the e 
sults as computed from them are within the range of accuracy required. aa 
_ Since exp Up T ¢ and exp Ln T become large for large values of Un Tf, 
the number one in the numerator and in in the denominator of Eq. 6 ‘may bo both — 


i 
| 
= KD (2)(2) (1 +b_) f exp-uy Td7.... 
| 
i 
q 


> 


For small values of | T, ap 


0 lim a a, * 


| Fig. ne of . shows curves of an versus Bn T for each of the constant val- 

= 0.2, 0.4, 0.6, and 0.8, - Inspection of this figure shows that, for all 

four ‘curves, ay ty/T for T 0.01, the lowest value for which an has 

. ye The Water Table for Low Values o fe ee —It may be s shown that the water _ 

Ho 

- table would not fluctuate in response to intermittent recharge if ug T— 0. 

imit problem (a hypothetical one from a physical point of view) reveals — 

 daanent features of the behavior of bad water table when 1. T is small, and 


the method of proof will be indicated. _ i 
The limit value of the ‘general term of the infinite series in Eq. 9( (a) is 


4 lim (2 n+ n+1)73 (-1)" exp - te? ((2n+1)2y_ 


\q 


it found that 


7 same result may be obtai 


charge | the recharge - which t1 is thus found, from 


> \ 


It is noted that for - 


t of the water table above the water level 
x in the drains for steady state flow ae. M.6), the continuous rate of recharge 
The derivation that, for small values of Ho the water 


table 1 will have i an approximately parabolic shape when the ground water is 


ac 
n 


35) 


- x), for 


| 
which N is a numerical value of pred 
— 
4 
' 
| 
| 

7 in wh | 
, 

q 
bo 
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R WATER TABLE. 
_ recharged intermittently and dynamic equilibrium has been reached. It also. 
follows that the amplitude of fluctuation will be small compared to the aver- 
age level of the water table, a feature brought out by numerical evaluation of | 
The transient Eq. 11 also leads to Eq. for up T From Eq. 2. is 


jen for the limit problem here under consideration that there would be no- B's 


rise of the water table. This result indicates that it takes many periods T to — 
‘approach equilibrium if 1. T is small. The latter conclusion has important — 
‘The previous considerations are not dependent on the value of ty. In fact, 
the factor ty does not occur in the final result (Eq. 37). The preceding con- 
clusions are therefore also valid for ty — 0, that is, for intermittent instan- 


Double Periodic Recharge. —In order to evaluate effect of irrigation 
deep percolation losses on the water table, it is convenient to divide each year 
- into two periods; the irrigation season and nonirrigation season. For the 


» irrigation season, it wilt ania that the ground water is recharged at 


~ 
! 
ie “actuat | HYPOTHETICAL RECHARGE (POSITIVE or NEGATIVE ) 


FIG. 2,THREE RECHARGES 3ES SPACED A ar INTERVALS OF T DAYS REPEATED 


regular intervals with every irrigation or each alternate 7 
q 


ing the nonirrigation season, there is presumably no recharge. 
‘The problem of ‘double periodicity’ requires an . extension of the theory of 
intermittent recharge, which is most easily explained by first analyzing the — 
example shown in Fig. 2. It requires superposition of Eqs. 2 and M.38. The © 
_ solution to the problem is obtained by superimposing hypothetical positive and 
“negative recharges as explained in ( MM), 2» and illustrated in Fig. 2. 
tion is drawn to an error in the first paragraph of the second column on p. 
3 557 of (MM). The second sentence reads “Let (r- - 1) recharges be applied 
which should be corrected to read “Let r - recharges be applied. ..”). 
In Fig. 2, three instantaneous recharges q are shown to be applied « each year, 
spaced at 28-day ( = T) intervals, so that, for practical purposes, , each year 
may be divided into thirteen 28- day intervals. Intermittent instantaneous re- 
charge is depicted in the figure, but both intermittent instantaneous recharge = 


“and constant — will be considered in the 


| 

7 
— 

a 


‘tal total year be expressed in a multiple of this interval. The latter is always 
possible because slight shortening or extension of the year does not signifi- ay 
affect the result. _ The general problem will be treated later. 
os _ The problem shown in Fig. 2 is defined by Eq. 21 for values of § < 3 (in a 
the first year). For 3 <  <13, the water table level is defined by Eq. M. 38. mig _ 
/- the second year, we need to evaluate the effect of the new recharges in \that 7 
year as well as the residual effect of the recharges in the first year. In the q 
= year, the effect of the recharges in that year 1 needs to be evaluated in i - 
~— eonjunetion with those of | the second and first years, and so forth for subse- 
que years. Water table rises will first be determined for the intervals 4 
 2+13(m - 1)< €<3+13(m- 1), i in which m = 1, 2, represents 


_ water table level in the first year, and for the preceding range of £, 


in which. the subscripts p, t, x have their - usual meanings; represents 
" the value to be used in the periodic function, and (2 +y) that to be used for 
“2 in the transient. a The subscript 1 of hy, 1 on the left hand side of Eq. 38 is 


is indicate that this water table level is ‘that for the first year of recharge. _ 
_ The effect on the water table in the second year of the three > recharges — 


— — the first year is computed from Eq. M. 38, that is, non 


hex 12 Dex 15+ ty ee ee .(39) 


- sidual height of the water table in the ‘second year, due to the recharges in 
the first year; the subscripts (12+y) and (15 +y) on the right hand side 
values of £ to be used in Eq. 11 or Eq. M.30 for the transient Se 

tI general, the effect < on n the water table level in later y gies of the three 
applied during the first st year may be 


wh 


for any year except the first year. ‘This occurs because 1i 
cludes the periodic function in which the variable y can only | vary over a unit a 


¥ interval of £, that is, from 0 to 1 (comment below Eq. 4). The height hy », 
‘. _ decreases monotonically, that is, its numerical value decreases. as m in- 
— 4 ‘creases and equals zero after an indefinitely large number of years. For the 
usual drainage situation, it is negligibly small after only afew years. _ 
| The resultant water table in any year, due to recharges in successive 
years, is obtained by simple Addition of the water table rises as computed 


from Eqs. 38, 39, and 40. For the mth year it follows that 
validity of this ‘equation may be ex) explained by redefining the symbols 


’. 9, and so on. The symbol hy 1 represents the water table rise due to re- 
charges in the last mth year, hy. 2 is the water | table rise due to ao in 


The “methodology to be developed require t, for each particular prob-— 
: 
— 
3 
| 
7 
— 
q 
3 
7 
W 
— 


WATER R TABLE 8 
due to recharges in the first maximum height of the wa water tutte: 
4 ~ the mth year is obtained from Eq. 41 by putting y = ¢ for intermittent constant - 
and by putting y= 0 for intermittent instantaneous recharge. 
_ Eq. 41 was derived for values of ¢ between 13(m - 1)+ 2 and13(m-1)+3 
a may be extended to cover each whole irrigation season, that is, values of E : 
"may vary from 13(m. - 1) to 13(m - 1) + 3. The variable y in the periodic. a 
= then goes from 0 to 1 as ¢ increases with each irrigation interval 2 
= one (positive) integral number to the next. : 
‘Similarly, : we find for the resultant water table h’ x,r dur 


The symbol ra a) 1 of Eq. ore the water table rise due to recharges ~ 
in the mth year, by, 2 is the water table rise due to recharges in the (m - i 


in EQ. 42. It may also be taken to indicate the — that is to be used in any il 7 
term (including the first term) of the series on the right hand side of Eq. 42. 
_ Any additional year requires one additional term in the series. The a P 
: height of the water table at the end of t the first nonirrigation season is obtained — _ 
: from hx, 1 by putting § = 13. The resultant water table (hx, r) at the end of the | 
 £=13(m- 1) for hy, (m- 1) , +++, £= (13) (2) for and, finally, =13 
_ ‘The preceding values of — may be generalized to any number er of recharges — 
plied and let the year be divided into b intervals of T days each. The irriga- n 
tion season, or more appropriately, the season of recharge, then covers - 
intervals, and the nonirrigation season, the season without recharge, covers 
a) intervals. Furthermore, let w’ the ‘time’ within the year expressed 


mth season is computed from 43 by using = 13m for hx. 
sg the irrigation season spaced at any interval. Leta recharges be ap- 


’ further into wt w represents the integral numbers 
_ from 0 to b, while y varies from 0 to 1. _ The latter subdivision prevents con- 
fusion as to which value is to be used for y in the periodic function (Eqs. 5 or 
M. 29) concurrently with the £’s for the other terms. For the irrigation sea- 


son and for any mth year but the first ‘year, the fol natin result is obtained: : 


+ hy (m-2)b- atwty ~ Mx, (m-2)b+wt+y 


7 
g 

| 

q 

| 

in unltS » lat varies Irom VU | ginning O e irriga 5Ca- — 

q 


the nonirrigation it follows that 


2)b- - tx (m- 


Observe that all (except the first) terms s of the right hand side of ne 44 
and 45 are identical in form. Furthermore, since the subscripts (b-a+w _ 
_.-® y), and so | on, represent values of E, it follows from Eqs. 11 and M.30 that 
the terms of both Eqs. 44 and 45 each form two geometric progressions (ex- 5 
cept the first term of Eq. 44 for the periodic function) which may be summed 
separately. Using Eqs. 11 or M.30 and performing the summations | 
Eq 44 yields, for the irrigation season, the following results: 


gle 


of which Nox 


= tate is obtained by putting w =a - 1, and 


1 


2n+1 
_/x exp Hy Ta - 1 
|exp a, Te =1| - Tw ty)... 
of which hpix is ‘computed directly from Eq. M.29 and w+y from Eq. 
_ M. 30; the maximum water table is obtained by putting w = -a-1, andy=0. 7 


"The factor bp, is given by Eq. 23. Eqs. 46 and 47 are valid for 0 < (w + da 


a. The ‘notation in in Eqs. | 29 and 0 should be changed as" indicated d by 
Eq. 


Similarly, t 
1. 


For nial it constant ant recharge (Eas. 5 5 and = 


3 


| sin (2n 

a(® Ja 


48 
(48), 


1961 m2 | 
: = = in which 0 S$ (w + y) < a. The maximum water table for any year is obtained [7m 
by, ting 4) = =, in ten onstan recharge or y = 0 
&§ 
| 
i | 
| 
11; the | 
on, from Eq. 49: 
| — 
_ 
= &g 
— 


WATER 
of which the first two terms are from 1. 11; the minimum 


- Uy T(w + y) 


_ Note that the notation on the left hand side of the preceding equations has 
son n changed from hx,r r and hy r r to hx, ¢ e and hx, e , respectively. ‘This is done > 
because Eqs. 46, 47, 48, and 49 describe an equilibrium condition (theoreti- > 
cally for m—+) equivalent to the so-called dynamic equilibrium defined by 
the periodic functions in (MM) for ‘single periodic’ recharge. . The transients" 
have been omitted ‘because the ultimate equilibrium water table level is of 
By interest in open ditch or tile drainage problems. The water table 
— levels “midway between the drains (he and he) are determined : from Eqs. 46, 
4, 48 and 49 by putting (- -1)" for sin (2n+1) 7 (x/L). Forw+y =a, Eqs. _ 
_ 46 and 48 are identical and so are Eqs. 47 and 49. Eq. 46 withhw+y=0is @ 
A identical to | Eq. 48 with w + “7% b. _ Eq. 47 with w + y = 0 is identical to Eq. 49° 
with w + y = b if the water table rise due to one recharge (q/e) is added to 7 
pay Equations for the rate of discharge are not given but may be obtained in © 
exactly the same manner. _ In fact, the equations for ‘double periodic’ inter- a 
_ mittent constant recharge may be immediately written down from Eqs. 14(a) 


18 from om inspection of | the foregoing — 


4 al 


Of all the equations presented in the paper, those for intermittent instan- 


a taneous recharge are the most useful for practical work. Their greater sim- 

 plicity makes them attractive. Furthermore, numerical evaluation shows that _ 

- there is often little difference between water table levels computed for an § 
interval (of recharge that is or zero, providing the amount of re- 


a irrigated agriculture, recharge to the groundwater approaches the condition J 
_ Of instantaneous recharge in many cases. The word ‘approach’ is used ad- 
visedly; instantaneous recharge, as defined mathematically, obviously y can 
never take place in nature. Field observations, for which data have been | 
_ gathered by Dumm, have | shown that this recharge condition is often closely 7 
approximated in normal irrigation "practice. Comprehensive > theoretical 
‘studies on flow of water through unsaturated soil support the field observations. — 
are, of ‘course, for which the of recharge must» 


drainage if one wishes ‘to. the height of the 


— 
4 
7 4 
q 
of which the fi 
mum water table is obtained by putting (w +y) =a. Eqs. and are valid 
4 
? 
— 
— 
| 
fluctuations 
| 


“1961 


‘oo conductivity, _ | The e effect on the water table of an overburden 
with a low hydraulic: conductivity in the vertical direction has been ignored - 
in the present study. The solution of that problem requires an extension all 
_ the theory to include evaluation of hydraulic head losses in the slowly pe perm- 
stratum. . It then to two simultaneous 


_ Surface 


“Mmpervious Borrier 


.—FLOW THROUGH A HORIZONTAL HIGHLY PERMEABLE LAYER 


Ane 


tial is done by to. determine22 the effect on the 


_— of the work ‘reported ir in this | paper was done by tl the w: writer while he 
- was employed by the U.S. Bureau 0 of Reclamation, Dept. of Interior (USBR). 
‘Earlier technical memorandums on on this ‘subject prepared by Glover were of 
great value to the writer in his study of the problem. The encouragement of 
‘Messrs. P. H. Berg, M. ASCE, and J. T. Maletic is gratefully acknowledged. 
Conversations Movers. Glover, Dumm, and Tapp were stimulating and 


_  22*The Tranmission of Tidal Waves in Elastic Artesian Basins,” by J. Wesseling, 
Netherlands Journal of Agricultural Science, Vol. 7, 1959, pp. 22-32. a eee 


— 
{ water table in a slowly permeable overburden (Fig. 3) rather than the piezo- : 
> 
= 
FIP? 
q 
— 


a The “symbols adopted for. use in 


eference and for the use of discussers: 

aos = total number of recharges during one recharge season in problems _ 
involving ‘double periodic’ recharge; | the — of the recharge sea- 


= number of intervals T into which each year is divided in problems © 


Verve 


= weighted mean depth of the profile through which the flow occurs; 


= 1/2qL, of released p per period T from one 
side of the drain, after dynamic equilibrium has been reached; 
hp, Np, Hpi,h, he, = height of the water table 
cert between drains, measured from the line connecting the water 
levels in the drains; the subscripts 1 and h are used to denote the - 
minimum and maximum water table 


hp, x hpi, x)= height | of the water for any x, as computed from the 
periodic functions for ‘single periodic’ ‘recharge; 


(or hej x)= height of the table any x, as computed from the 


transient functions for ‘ periodic’ recharge; 


(or hy x)= height of the water table for r any x, measured from the line | 
connecting the the water table levels: in the edrains; a. 


-—— be x)=: = resultant height of the water table after a number of years — 
of ‘double periodic’ recharge; the symbol without prime denotes the 7 
water table during the recharge season, with prime it denotes the water 
table during the non-recharge season; 
_ (or hy m) = the rise of the water table in the mth year due to recharges 
in first year, as computed for ‘double periodic’ recharge; the symbol 7 
_ without prime denotes the water table rise during the recharge sinned 
prime it denotes the rise = recharge season; 
(or h x hix. height of t the water table as computed for the = 


condition (of ‘dynamic equilibrium from the periodic fu 
‘double p eriodic’ recharge; 


= constant rate a surface- recharge p per unit of land 


‘surface; 
= ‘p/KD, recharge coefficient for L D> hm and K K; 


width of aquifer or ‘spacing g of drains; 


i 
: 
| 
= 
| 
= 
7 | 
E 
a | 
|. 
£ 
mt 
— 


ty length of drain computed for flow nto th the drain from only one side of 


A 
— 


Tn = volume o ot groundwater stored below the water table a 


taneous recharge; 


interval of recharg within ‘the eriod T: 


= horizontal coordinate; a 

drainable porosity or or 


(K ‘D/e) (2n +1)? 


14159; 


time within period T; =: 7 <T; and 


Ho T= (K D/e) (1/L)? T, the response coefficient 


w + ‘time within ‘the: year “expressed in units of T; w’ varies from 


of the to (at end o of non-recharge 
‘The symbols y, £, are defined | by Eq. 2, aly) ply) by Eq. 3, an by 
. 5, and by, by Eq. 21; they are all dimensionless. The subscript i is used 1 


instantaneous ‘recharg eis assumed. 


= 
= 
stan- 
a 
= 
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WORLD PRACTICES | IN WATER MEASUREMENT AT TURNOUTS?@ 


= 
Closure by Charles W. Thomas _ 


and measurement of water at the farm turnout. 


EVALUATION TABLE FOR 


Operating | Relative Type Head Ran e of Expected 


nance ‘lequigment 


= It is regrettable that time and space die not "permit writer to evaluate 
each of the in the paper. A very useful contribution to 


aves, information as: basic. operating whether used in con- 
_ junction with, or independent of a control device, range of flow rates that can 
a. be measured, head loss, cost, and so forth, should be very useful if used with © F 
care. The headings for such a table might be similar to that given in Table 1. 
Perhaps this suggestion could be the stimulus for someone to expend time an Fs 
effort to produce such an 1 evaluation or charts. 


ulation, and measurement of flows. As mentioned in the paper a shutoff is, = 
most instances, ‘required at a farm n delivery because of the manner in which 7 


a June 1960, by | Charles W. Thomas (Proc. Paper 2530). ss 


18 ‘Hydr. Engr. ., Bur. . of Reclamation, U.S. Dept. of the interior, Denver, Colo. 


CHARLES W. THOMAS, 18 F. ASCE.—The discussion by Chow serves toil-_ 
| 
— 
| 
4 
a (l) | “| 
— 
ox! 


Many 
. it r may be e economical o or convenient to utilize control devices and structures as" 
indications of the quantity of water flowing. Such structures are for the most 
part designed without thought of their use for t this purpose. . The difficulties of 


be adequate for eqeentions! purposes and, in some instances, may be of a high 
order of accuracy. It must first be determined how well the control will — 
_ measurement and the } expected accuracy of results. A calibration may be 


obtained by application of data derived from similar devices, by hydraulic mod- 2 
el tests. or by measurements made inthe field. _ Regardless of the method a 
*% to: obtain a calibration of a control | device, the > accuracy of the calibration can 
a no better than the accuracy of the device or means used to effect the cali- - 
bration. It may be seen then that caution is indicated in making such nome 
tions, and and the results should be critically judged before general busage. 
_ ny operational Lape can be removed if a single device can be made ~ 
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systems are operated. Economy in installation can be effected if the 
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_ WILLIAM Ww. DONNAN, 2 F. ASCE. — Of particular interest to this writer 
Z to | follow the progress _and note that only by group action was it ‘possible 


to plan | and finance the engineering works required to bring water to the ia 


_ Angeles, California,area, 
5 


_ California laws and policies create a favorable climate for the ied 
ity of water districts to cope with and resolve water problems. ‘Today (1961) there 
. are over 2,500 municipal, commercial, mutual or other organized water os 
ice agencies in California and over 1,000 are located in Southern California.3 
ial If ‘Southern California is to ‘continue to expand in population and near 
growth, additional water must be developed and imported. The State Depart- _ 
_ ment of Water Resources estimates that, at the present time (1961), , there is - 
a a need for about 400,000 acre ft of additional w water per yr. 3 In looking to the 7 
_ future, if all the agricultural, urban and industrial expansion were to ean 
to an ultimate level, the mean seasonal requirement of ae water would 
be 4,000,000 acre ft. 3 How can this challenge be met? 
= _ The solution tothe water problems of the future rest inthe wikia of South-- 
ern California’s people to cooperate. There are a number of different schemes 
being investigated to developnew sources of water. Eachonehas certain mer- © 
its but each one will require the —e in cooperation, grouping of effort and © 
tat The conversion of sea aie involves the purification of sea water so 
that it canbe used for domestic and irrigation purposes. The best known meth- 
ods utilize tremendous amounts of electric current. - Moreover, any water de- 
velopment | will have tobe pumped from sea level to its point of use. At present 
z only large- -scale, high initial investment plants offer the promise of aneconom a 


. the Feather River. The proposed capacity of this reservoir would be 3,500, 000 
ft and the dam would” be 710 ft high. There would be installed electric 


the Sacramento River and then be pumped up to canals and conduits flowing 
southw ard toward Los Angeles. ‘The entire scheme would cost about 1.5 — 


& December 1960, by Samuel B. Morris 
B 2 Agric. Engr., Soil and Water Conservation Research Div. , Agric. F Research Service, © 
3 “Water Utilization and Requirements of California,” Bulletin No .. Publication of _ 
the Calif. State Water Resources Bd., June, 1955. 
0 “The Feather River Plan, 1,” by A. D. ‘Edmonston. | Presented in March, 1956 to the 4 
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dollars. It would require that Southern California cooperate fe with Central Val- : 

ley and San Joaquin Valley people to vote bonds for this great enterprize. The 
stages of the Feather River Project are started. 

Fa 3. The California Water Plan5 involves the utilization of the total soil and 


“. _ water resources of the State. The plan calls for building 260 new major res- q 


-ervoirs andan aqueduct systemto transport the water fromthe northern areas 
of excess to the arid south. It would cost an estimated 12 billion to 13 billion 
dollars. The Feather River Project would be only one facet of the total plan. 4 
This plan would require the utmost in cooperation between | areas within the 
State. For example, northern water rights would have to be safeguarded, more | 
metropolitan water might have to be sent to San Diego County and replaced with 
Feather River water, and water rights would have to be traded or or 
On November 8, 1960 the people of California voted bonds. for ce develop-_ 
“<i ment of the California Water Plan. . Fortunately, California has an outstanding — 
State Department of Water Resources. _ This: technical staff can be counted on 
a to ) provide all the skill and engineering science to develop these schemes. __ 
_ The big question is, can the people, all the people, organize and cooperate 
-& to bring these schemes to fulfillment? History has shown that only when two 


or more farmers join together can they divert a stream; only when a group 
organizes can they | build a reservoir; only when a large city cooperates can 
; they build an aqueduct. _ Only when cities and groups of cities combine can they 

build a Colorado River Aqueduct. 


Se of these problems might be facilitated by the creation of a water 


sources are e taxed for the development of water. This agency, working 
<< the State Department of Water Resources, could negotiate with and join with 
Northern ‘California groups to ultimately develop the total water Supplies 


CHARLES H. LEE, 6 F. ASCE.—The author presented, in a most inter- 
esting manner, , the early history of the water supply of the City sof Los Angeles. 

_ This history is of particular interest to the writer as he was, for many years: 

— (1906 to 1935), intimately connected with the activities of the Los Angeles Wa- 
ter Department especially with the Owens River Aqueduct. 
Los Angeles, being in a semi-arid region, has had to seek water r from great 

_ distances to supply the needs of its rapidly growing population. As stated by — 

the author, the early policy of the Water Department was not to develop nearby © 

_ sources of water needed for the growth of neighboring communities. _ This led 4 
tothe selection of the Owens River, at a distance of 250 miles, ‘and. later the 
_ Mono Basin, at a distance of 338 miles, and the Colorado River, 350 miles a- 


consideration of the use of water from the Feather River in Northern Califor- 
nia, a distance of 420 miles. It is noteworthy, however, that this policy has 


; a way,as sources : of supply. Continuing increase in demand is presently causing © 
been extended to within the of the Fernando Val- 


pioneered many new ideas in 1 engineering research and | {prone wig as well as” 

5 “Preview of the California Water Plan,” publication of the State Water Resources 

6 Senior Partner, Lee and Praszker, Cons . Civ. San Francisco, Calif. 
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construction the Owens River Aqueduct period, the use 


 -tufa cement was developed and practically applied. Much was also contributed _ 
to the mechanics of powered transportation over rough terrain. Early in 1908, 
Chief Engineer William Mulholland d initiated the first quantitative ‘survey of in- : 


ley where he recognized the great potentiality of ‘the Independence Basin as 
reserve source of supply for the Aqueduct in dry years. The writer was priv- 
ileged to conduct this survey, that extended over a period of several years, 
and led to extensive well development, both in the Independence region and _ 
upper portion of Owens Valley in the Big Pine and Bishop regions. During the > 
dry period between 1920 and 1930, this source was of material aid in main- 7 
taining the Aqueduct ‘supply. The basic hydrologic methods followed in making 


al 


this survey have today been developed into standard, more or less routine, — 
methods of arriving at hydrologic inventories and safe yield values for closed § 
iu q ground water basins. Such methods are in use not only among private engineers 
but also the engineer-geologists of the Ground Water Division of the U. S. Ge- F 
ological Survey and of various state agencies, such as the Water Resources | 
Department and the Water Rights Board of the State of California. ‘ae 
_ Another field of hydrologic engineering pioneered by the Department was 
the making of snow surveys in mountain watersheds for the purpose of estima- _ 
ting available seasonal water supply. This work was initiated about 1911 by « 
J. E. Jones, M. ASCE, Chief Hydrographer of the Department of Water and — 
- Power, and technique was developed in cooperation with other agencies working 
on the east slope of the Sierra Nevada Mountains. Such surveys later became — 
{ routine. practice throughout the Sierra Nevada and other high mountain ranges , 
of the West. In California and other western states, snow surveys are now 
conducted annually by State ite agencies, supplemented in ‘some areas by the Fed-- ; 


mot aids that have placed it among the leaders in the waterworks | field, 
The paper ends with the statement that the irrigation chapter within the City 
of Los Angeles is about to close. The impact of this statement on the ground 
water supply of the San Fernando Valley is notable, ‘because urbanization with | 
Bar areas and roofs is taking the place of irrigation. Furthermore, demand 
£ arisen for flood control on streams traversing the valley. These improve- 
ments greatly reduce the > opportunity for natural | percolation tothe water table 


= runoff in stream channels that have now been largely sealed by impervi- | 2 
ous concrete lining. . The effect of these aneimeees in depletion of the nat- 
safe yield of the basin is substantial. 

The author is to be congratulated for his neaeeahtee recording of the his- -— 


torical highlights of one of the outstanding waterworks systems of th the ‘United 


‘oe pe,” ” by Charles H. Lee, Transactions, ASCE, Vol. LXXVIII, 1915. ‘Sila ti 
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Discussion by George H. 


GEORGE H. HARGREAVES,!5 F, ASCE. — —Many methods are inuse for com- 
puting evapotranspiration from climatic data. Most of these, although useful a 
7 as in certain geographic areas, are either approximations or limited in suitability 
; . iG general application. Three methods appear to offer the advantages of ac- 7 
curacy and general applicability. These are the Penman Method, the black bulb 


pan offers the ‘advantages of 1 more e general availability of data and direct rela- 
_tionshipto evapotranspiration. A summary of some of the literature describing 

_ this relationship is given by Pruitt12 in his introduction. — 
Available research data indicate that irrigation requirements or consump- 
tive use requirements depend ‘primarily on the degree of ground cover and on 

_ the rate and stage of growth. Rapidly growing crops are e usually dark green 
and consequently are able to absorb a maximum of solar energy. A dark green 
7 _ crop providing complete ground cover will have a maximum potential consump- | 
ei tive use or potential evapotranspiration. This concept of potential evapotrans- _ 
piration is used by Penman, Thornthwaite, and Pruitt. 

order to relate climatic factors to consumptive use or irrigation r require- 
-ments, it would seem desirable that the | correlation between the factor and po- | 
tential evapotranspiration be evaluated. ‘Pruitt 1 measured potential evapotrans- 

- piration by clover at Prosser, Washington and correlated the results with evap- 
oration from six types of pans. Coefficients of correlation for the evaporation 
pans varied from 0.972 to 0.990. The coefficient of f correlation wit with black bulb 

s _ atmometer evaporation was found to be 0.972 and with Penman’s formula (Eo. 
- X 0.97) to be 0.986. For the Blaney-Criddle method (K = 1.08) the coefficient — 
= found to be 0.908, and for Thornthwaite’s method (E T x 1.72) it was found | 
be 0.849. data and discussion presented by Pruitt clearly support the 
use of pan evaporation data in determining consumptive use or irrigation 1 re- 

ce it is concluded that pan evaporation data give an naccurate i index of the | con- 

sumptive use of water by agricultural crops, then a comparison of pan evapo- 
ration with results from any given method of computation sien a means of : 

evaluating that particular method. 
Munson’ s method of estimating consumptive use temperature data by 
use of Table 1 can be evaluated by using a climatic atlas in order to obtain - : 
= visual representation of climatic conditions. Fig. 2 shows temperature 
a and evaporation conditions for the United States during January. Fig. 3 shows © 
a December 1960, by Wendell C. ‘meen (Proc. Pap 
15 Irrig. Advisor ICA, Rio de Janeiro, Brazil, 
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conditions during July. The climatic data —_— inthe figures are from Stephen _ 
S. Visher.16 Evaporation data used in Figs. 2 and 3 are from reservoirs and 

‘Shallow lakes and are approximately proportional to pan evaporation. Fig. 4 _ 


shows normal May through October evaporation from pans and normal summer 


_ Blaney10 ; shows monthly ratios of consumptive use of water by alfalfa to pan 
4 wimeminn. These vary from 0.45 to 0.86 with an average value of 0.66. Low 


rates of use by alfalfa occurred after cutting of the crop thus reducing thetotal | 

iM water requirement . Pruitt!2 measured the consumptive use of water by Ladino 
clover on two plots cut at staggered times so that one plot would always pro- .- 

100% ground clover and consequently full potential evapotranspiration. 

a _ Consumptive use or evapotranspiration under these conditions was found to be 

a 7 P 108% of evaporation from a U. S. Weather Bureau Class A pan. By selecting a 7 


_ suitable factor to represent the degree of ground cover, evaporation, as shown 
in Fig. 4, provides a convincent index of the consumptive use of water for ag- I 


__-Errata.—Two typographical errors should be corrected. On page 54, line 6 6, 7 
black atmometers absorb 
7 At the a of Table 7 (page 55) the reference should be 12 instead of 10. 


“Figs. and 4 indicate that variations in evaporation of of 
consumptive use are considerably more complex than variations in tempera: 
a ture. An inspection of Fig. 3, indicates a fairly large area inthesouthernpor- _ 7 
of the United States with approximately the same mean monthly tempera- 
ture. Evaporation within this area of approximately the same temperature var-_ 7 
ie ied trom a low of 5 in. in the more humid areas to a high of 12 in. or more in i : 
the arid areas. Although temperature may provide a convenient and useful in- 
a dex of the consumptive use of water for some areas, , the method proposed by _ : 7 
‘ 
Munson should be used with considerable caution in others. 


more than 90% of the solar energy instead of 50%. 


-@ sree a 


am “Climatic Atlas of the United States,” by Stephen S. Visher, Harvard Univ. Press, 
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Discussion | by Vahe Sevian and Kar Kamil 1 Taj- -Eddin 


_ VAHE J. SEVIAN, 3 F. ASCE and KAMIL TAJ- EDDIN. 4_ Varieties of wheat 

and barley 1 were cultivated in the northern regions of Mesopotamia 60, centu-— « 
ries B.C. Agricultural life spread in central and southern Mesopotamia 20 | 
7 centuries later, whereas after a further period of 20 centuries (bringing us sd 2) 
cirea 2000 B.C. ) greater varieties of land products were obtained. Rice eel 


or the “Hanging Gardens of Babylonia,” ot prosperity of the valleys continue 
and attained its height under the brillant achievements of the Abbasids (754 
= through 1258 A.D.), , ending with the great blow given by Hulagu in 1258 A. Dd. 


At the turn of the present century, efforts were noticeable for improving ir- 
agriculture. Progress, though slow, was appreciable. From mid- 
century onwards greater stress is made in improving the agricultural pattern 


7 through modern technics, by making better use of the water, by restoring to — 


the soil in many regions its lost productivity, and without neglecting the social 

_ aspects of the various problems involved. 
The: author mentions that “successful development of irrigation in the val- 

- ley... is greatly dependent on good engineering, progressive laws... and y 


~ cooperation with neighbouring countries. ” These are obvious; an important 
point is the need of an integrated program of development of the main re- 


sources: land and water involving human association. 


Development of water resources will imply with countries 


on which the valley depend for their supplies. These were not needed during ce, 
past centuries, for the basin was almost entirely within one country, the Otto- : 


m prevailing then. i With the 2 end of First World War, Syria and Iraq ceased to be 
parts of the Ottoman Empire. Consequently, the necessity for understandings — : 
is now felt with serious considerations given for agricultural and other develop- 


Since July 1558, the country is amidst an “Agrarian Reform” process. This 
_ will change the thousands of years practices and land d (agricultural) tenure in- 

to a new aspect more adaptable to advanced conditions. 
we Throughout history, the “Great Flood” was followed by very many others. _ 


Their destructive effects were byt the extent of done to. rural 


@ December 1960, by Banley S. Butler (Proc. Paper 
3 Formerly Inspector General of Irrig., Ministry of Agric., Baghdad, Iraq. 
4 Engr.-Secry. to the Dir. Gen. of Irrig., Baghdad, Iraq. —™ 
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and urban settlements in ntheir varying of advancement. recent. years 
_ the Euphrates, in 1940, and the Tigris, in 1946 and 1954, caused very serious — 
_ flood damages estimated at $200,000,000 for the 1954 flood mentioned in the 
1941, floods on the Euphrates valley are 2 sufficiently ‘controlled by 
works at Habbaniya. . On the Tigris flood risks are almost eliminated ted through — 
Tharthar project. completed just after the 1954 flood. 
Tidal irrigation has for centuries been practiced inthe estuary. of the rivers. 
“ad Further upstream, gravity flood and flood irrigation was applied. This process 
is still practiced but greatly improved by means of water diversional structures | 
- and methods. There is another method of irrigation, with over 5,000 small 
4 individual pumps erected along the rivers and some canals. In the northern © 
semi- humid regions, irrigation depends on the vagaries of rainfall. > 
2 A method timidly used till now is the utilization on a greater scale of the 
ean water resources wherever suitable in quantity and quality in regions — 
distant from main streams. Investigations and surveys were undertaken and — 
the author’s contribution in this respect is acknowledged. 
As regards to ultimate possibilities of irrigation with reference to surface 


and ground waters, it seems advisable to ascertain first volumes to be with- 


drawn fromthe rivers outside Iraq and then to decide onthe future agricultural 
_ pattern and probable change either totally or regionally from the present fal- 7 
low to intensive cultivation, 
Dam sites exist in Iraq and in the headwaters of the rivers. rs. Though the 
_ depression of Abu Dibbes on the Euphrates and Tharthar on the Tigris have 
been mentioned for flood control and storage, it will after proper investigations 7 
appear that other possibilities will afford better results for water conservation 
On the Euphrates, Turkey’s demands will mostly be for hydro- electricity. 
‘The retention of the water in the reservoirs and subsequent ‘releases may not 
suit the agricultural requirements of distant Iraq, unless such releases are 
harmonized through a reservoir preferably in Iraq. In Syria irrigation re- 
‘quirements have priority over hydro- ‘electricity to a certain extent. 
wo The author mentions - from personal knowledge - - that “in several areas a 
bordering Iran there is often a severe shortage of surface water for irrigation | ' 
y in Iraq.” The water is used for areas across the border with disregard to the 
a interests and vested rights from time immemorial of the frontier regions of 
_ Iraq. Understanding between the interested countries is essential, ground wa- 7 
ter will not satisfy the requirements of the affected areas, only water diver- 
sional works from the Diyala river will help in satisfying the region. Rect 
é the flow of a stream passing through Khanagin (a frontier rural and urban cen- 4 
tre) was entirely cut before crossing the frontier 
7 tributed in canals controlled by the Irrigation Department in strict proportion _ 
3 to cultivable areas. In some ancient canals still in operation, old established 


rights are maintained for the time being. Years are to — a 


Pe in many fields has been asiheved in a small region of the world — irrigation 
_ has been practiced for over 60 centuries. Ce 


7 
&§ 
=| 
| 
Al 

| 

| 

| 
— 


Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
‘@raulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (BA), Soil 
_ Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterwrys and Harbors 
(WW), divisions, Papers sponsored by the Department of Conditions of Practice are identified by the sym- 
- bols (PP). For titles and order coupons, refer to the appropriate issue of “Civil Engineering.” 
with Volume 82 (January 1956) papers were published in Journals of the various Technical Divisions. To 
locate papers in the Journals, the symbols after the paper number are followed by a numeral designating 
the issue of a particular Journal in which the paper appeared. For example, Paper 2703 is identified 
2703(ST1) which indicates that the paper is 
JUNE: 2494(IR2), 2490(S 2497(EM3), 2501(SM3), 2502 
(EM), 2503(PO3), 2504(WW2), 2505(EM3), 2506(HY6), 2507(WWw2), 2508(PO3), 2500(ST6), 2510(EM3), 2511 
-(EM3), 2512(8T6), 25139(HW2), 2514(HY6), 2515(PO3), 2516(EM3), 2517(WW2), 2518(WW2), 2519(EM3), 2520 
2521(HY6), 2522(SM3), 2523(ST6), 2524(HY6), 2525(HY6), 2526(HY6), 2527(IR2), 2528(ST6), 2529 
(HW2), 2530(0R2), 2531(HY6), 2532(EM3)°, 2533 2534(WWw2), 2535(HY6)C, 258% P03)¢ 
JULY: 2542(ST7), 2548A4), 2 2545(ST7), 2546(HY7), 2548(8U2), 2540(SA4), 
9550(SU2), 2551(HY7), 2552(ST7), 2553(SU2), 2555(ST7), 2556(SA4), 2557(SA4), 2558(6A4), 
AUGUST: -2565(E M4), 2566(ST8), 2567(EMé), 2568(PO4), 2569(po4), 2570(HYS), 2571(E M4) 
2572(EM4), 2573( Bnd), 2574(sha4), 2575(EM4), 2576(EM4), 2577(HY8), 2576(E M4), 2579( POS), 2880 
EM4), 2581(ST8), 2882(ST8), 2583(EM4)¢, 2584(PO4)¢, 2585(ST8)¢, 2596(sm4)°, 2567(Hy8)°, 
EMBER: 2588(IR3), 2591(IR3), 2592(HW3), 2563(IR3), 2594(1R3), 2595(IR3), 


a, KWW3), 2607%(5A5), 2608(8TD), 260 2611(WWS)°, 2612(8T9)¢, 261 


OCTOBER: 2615(EM5), 2617(8T10), 2618(SMB), 2621(ST10), 


 2622(EM5),  2623(SM5), 2624(EMS), "2625(SM5), 2626(SM5), 2627(EMS), 2628(EMS), 2629(ST10), 
2630(ST10), 2631(POS)©, 2632(EMS)¢, 2633(ST10), 2684(ST10), 2685(ST10)°, 2636(SBA5. 
OVEMBER: 2637(ST11), 2638(ST1 1), 3639(CO3), 2640(ST 11), 2641(SA6), 2642(WW4), 2643(ST11), 2644(HY9), 

264608 2647(Ww4), 2648(WWw4) 2649(WW4),2650(ST1i), 2651(CO3), 2652(HY9), 2653(HY9), 
2654(ST11), 2655(HY9), 2656(HY9), 2657(SA6), 2658(WW4), 2659(WW4)°, 2660(SA6), 2661(COS). 2662(C08), 
ECEMBER: -2668(ST12), 2669(IR4), 2670(SMG), 2671(IR4), 2672(IR4), 2673(IR4), 2674(ST12), 2675 
(EM6), 2676 (IR4), 2677(HW4), 2678(ST12), 2679(EM6), 2680(ST12), 2681 (Sms), 2682(IR4), 
2684(SM6), 2685(ERA), 2606(EMS), 2687(EM6), 2688(EMS), 2689(EMS), 2690(ER 

ST12), 2603 (ST12), 2695(IR4)°, 2696 (sM6)°, 

“(HYL), 2707HY1), 2708(HY1), 270@(PO1), 2710(HY1), 2711(HY1), 2712(ST1), 2713(HY1), 2714(PO1), 
(ST1), 2716(HY1), 271 2718(8Al), 2719(SU1)®, 2720(SA1)°, 2721(ST1), 2722(PP1)°, '2733(P01)°. 2724 
BRUARY: 2726(WW1), ‘2727(6M1), 2728(EM1), 2720(WW1), 2730(WW1), 273 1 (EM1), 2732(SM1), 2733 (WW1), 

2735(EM1), 2736(EM1), 2737(PL1), 2738(PL1), 2730(PL1), 2740(PL1), 2741 EMI), 2742872), 
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(Hwi)°, 2783(SA2)°,.2784(CO1), 2785(CO1)®, 2786(IR1)©, 2787(ST3)°, 2788(AT1)°, 2789(HW1). 
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